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ABSTRACT 

We use the Marcario Low Resolution Spectrograph (LRS) at the Hobby-Eberly- Telescope (HET) to 
study the kinematics of pseudobulges and classical bulges in the nearby universe. We present major- 
axis rotational velocities, velocity dispersions, and /13 and /i4 moments derived from high-resolution 
(o'inst ~ 39 kms~^) spectra for 45 SO to Sc galaxies; for 27 of the galaxies we also present minor axis 
data. We combine our kinematics with bulge-to-disk decompositions. We demonstrate for the first 
time that purely kinematic diagnostics of the bulge dichotomy agree systematically with those based 
on Sersic index. Low Sersic index bulges have both increased rotational support (higher v/a values) 
and on average lower central velocity dispersions. Furthermore, we confirm that the same correlation 
also holds when visual morphologies are used to diagnose bulge type. The previously noted trend of 
photometrically flattened bulges to have shallower velocity dispersion profiles turns to be significant 
and systematic if the Sersic index is used to distinguish between pseudobulges and classical bulges. 
The correlation between /13 and v/a observed in elliptical galaxies is also observed in intermediate 
type galaxies, irrespective of bulge type. Finally, we present evidence for formerly undetected counter 
rotation in the two systems NGC 3945 and NGC 4736. 

*Based on observations obtained with the Hobby-Eberly Telescope, which is a joint project of the University of Texas at 
Austin, the Pennsylvania State University, Stanford University, Ludwig-Maximilians-Universitat Mnchen, and Georg-August- 
Universitot Gbttingen. 

Subject headings: galaxies: bulges — galaxies: evolution — galaxies: formation — galaxies: structure 
— galaxies: dynamics 



1. INTRODUCTION 



There is ample observational evidence that bulges in 
early type spiral galaxies come in different varieties. 
They are not all just like small elliptical gal axies which 
happ en to live in the centre o f a sp iral disk (jKormendy 
19931 [Kormendy fc Kennicutt||2004[ ) 
While classical bulges seem to lie 



et al.|2010 1 , although it seems ubiquitous in this class of 
objects. 

Internal secular evolution is commonly seen as an im- 
portant ch annel for the formation of central disk-like 
struc tures (Kormendy & Kennicutt 2004 Athanassoula 
20051. In this picture, asymmetries such as spiral struc 



on photometric pro- 
ject ions of the the fundam ental plane of elliptical galax- 
ies (Fisher & Drory 2010) pseudobulges resemble disks 
more than little ellipticals. They are still photomet- 
rically distinct from the outer disk as they appear as 
a central brightening above the inwards extrapolation 
of the outer exponential disk profile. However, as op- 
posed to classical bulges, their Sersic in dices fall close- 
to or below two (Fisher & Drory 2008). Other groups 
have shown that a large traction of galaxies with boxy 
or peanut shaped bulges d o show signs of inner disks 



ture and bars support the angular momentum transfer of 
disk material and thereby the transport of gas into the 
inner bulge regions. In t heir recent study of bu lges within 
the local 11 Mpc volume Fisher & Drory ( 2011 ) show that 
a majority of bulges in the local universe are pseudo- 
bulges. Their existence in large quantities in our local 
(low density) environment may seem to pose a problem 
for the un derstanding of the baryo nic physics of galaxy 



formation (Kormendy et al 
is not clear how 



.20101) as, at first sight, it 

central disks would survive the large 



(gurcau & Freeman 1999; Chung & Bureau] |2004 
[rmendy fc Barentine||2010| ) 
disks is however not limitec 



Ko- 



T'he phenomenon of inner 



quant ity of low redshift (z < «1) m ergers ( Stewar t et al. 
2008 |) typical of ACDM simulations ( [White fc Rees|1978 



White & Irenk 



imited to bulges that morphologi 



1991). But disk structures do not have 



cally resemble disks as a whole ( Scorza & B 


ender"1995 
2003, 2004 


Emsellem et al.| 2004| ll'alcon-Barroso et al. 


Sarzi et al.||2006| |Falc6n-Barroso et al.||2006| |Comer6n 
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to be destroyed m all mergers. [Hopkins et al.| ( [2008 1 
show that the heating of the disks m a minor merging 
event is a non-linear function of progenitor mass ratio 
once the satellite rigidity and the orbits are modelled 
properly. In addition to the impact parameters and the 
mass fraction, the baryonic gas content within the pro- 
genitors is an impo rtant parameter to the final result of 



a merging process (Hopkins et al. 2009 Stewart et al. 



2009[ 

likely to produce disks 



Governato et al. 2009[ ) — wet mergers are more 
Minor mergers may a lso create 



inner disk structures ( Eliche- Moral et ar[[2011 ) while in- 
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creasing th e Sersic index only moderately ( Eliche-Moral 
et al. 112006 1 . Based on semi-analytical models tor hierar- 



chical gr owth which include pre scriptions for the survival 
of disks, Fontanot et al. (2011) show that the existence 
of the majority ot the galaxies with no significant bulge 
component in the local volume can be understood in the 
context of ACDM. 



Kormendy & Kennicutt (20041 identify a number of 
criteria to ditterentiate between classical and pseudo- 
bulges such as the bulge vs. the disk ellipticity, their 
locat ion in the VmaxJc r diagram an d the bulge morphol- 



ogy. Fisher fc Drory (2008 20101 show that the bulge 



Sersic Index is successful in ditierentiating bulge types 
— pseudobulges have Sersic indices that fall near or be- 
low n = 2, unlike classical bulges and elliptical galaxies 
which have higher Sersic indices. 

The identification of the bulge morphology as well as 
the accurate derivation of photometric structural param- 
eters heavily relies on high spatial resolution imaging. 
Recently, such data has become available for a large 
number of bulges. Adding sufficiently high-resolution 
spectroscopic data allows us to ask the question whether 
all morphologically disk-like bulges also show kinematic 
disk-like behaviour, such as high v / a values and/or flat- 
tening or drops in the a profile? Also, whether differ- 
ences seen in structural parameters such as Sersic index 
are reflected in the kinematic structure as well? 

In [|2]we describe the sample selection and characteris- 
tics, in Sj3] we describe the long-slit observations, in S|4] 
we give account on the derivation of i/-band surface 
brightness profiles and their decomposition, as well as 
the details of the data reduction, especially the removal 
of emission features, and finally the kinematic extraction. 
Our results are presented in S|5) in S|6]we discuss correla- 
tions between kinematic parameters and morphological 
parameters and indications for an increased rotational 
support of pseudobulges. We finally discuss and summa- 
rize our findings in SJtI 

2. SAMPLE 

As we aim to study the kinematics of bulges, our sam- 
ple consists of 45 galaxies spanning the full range of 
Hubble types that do contain bulges: SO to Sbc. Fur- 
ther, roughly two-thirds of our galaxies are barred, a 
similar fraction to the total fraction of bars observed in 
the local universe (see Fig. [T]). Table [T] lists the ob- 
jects in our sample. For signal-to-noise (S'/A^) reasons 
we are biased towards high-luminosity objects. Absolute 
B-Band magnitudes span the range from Mb = —17.3 
to Mb = —21.3. Central velocity dispersions lie be- 
tween 60 kms^^ and 220 kms~^. We select our targets 
to be located close enough in distance to properly re- 
solve the bulge regions in typical seeing conditions. With 
the exception of NGC2964 and NGC4030, all galaxies 
have bulge radii larger than 5 arcseconds and are typ- 
ically located at distances closer than 25 Mpc. Only 
NGC4030, NGC4260, and NGC4772 are located at sig- 
nificantly larger distances of 29.3 Mpc, 48.4 Mpc and 
40.9 Mpc, respectively. The bulge radu of 7.3 arcsec- 
onds and 23.5 arcseconds of the latter two leave us confi- 
dent that we are able to nevertheless sufficiently resolve 
their bulges. NGC 2964 and NGC 4030 have bulge radu 
of 3.1 arcseconds and 3.0 arcseconds, and are excluded 
from all structural analysis concerning the bulges, we re- 



strict ourselves to presenting their kinematic data. 

In order to break the known degeneracy between the 
bulge effective radius and Sersic index in ID sur face 
brightness decompositions (Graham & CoUess 1997) we 
require all our targets to have MS'l' imagmg iii I' loOW, 
I, or R band. Most of the objects are found in |Fisher| 
fc Droryl ( |2008[ ) and/or [Fisher &: Dr"OTy| ( |2010[ ) and have 
extensive HS'l and ground- based multi- wavelength cov- 
erage. To allow for a visual inspection and morpholog- 
ical classification of the bulge region we select objects 
which have clos e-to y-band HST images of their bulge 
region (see \ 4.3 ) available from the archive and we do not 
observe edge-on or close-to edge-on objects (i > 70°). 
Four objects in our sample do have a larger inclina- 
tion. NGC 1023 and NGC 4371 are SO-types and con- 
tain very little dust and allow an undisturbed view into 
the bulge region. The situation is different for NGC 3593 
and NGC 7331, where the inclination and — in the case 
of NGC 3593 — the absence of an optical HST image 
inhibits a morphological classification. We present the 
kinematic data for those objects refrain from classifying 
them as classical or pseudobulges. 

3. OBSERVATIONS 

We obtain major axis spectra for all and minor axis 
spectra for about half of the galaxies in our sample. In 
some cases the observed position a ngle is not identica l 
to the one published in Hyperledapi] ( [Paturel et al.|2003 ). 
Also, in a few cases the minor axis position angle is not 
orthogonal to the major axis position angle. We label 
observations accordingly in Tab. [2] 

Observations were carried out in service mode during 
the period from April 2005 to Aprfl 2010 (see TahM at 
the Hobby-Eberly Telescope (HET; [Ramsey et aL|| ra98l 



at McDonald Observatory. W e use the Marca rio Low 



Resolution Spectrograph (LRS; [Hill et al. 1998) with a 
one arcsecond wide and 3.5 arcmmute long-sht, the E2 
phase volume holographic GRISM, and a Ford Aerospace 
3072x1024 15 ^m pixel (usable range 2750x900 pixel) 
CCD detector yielding a spatial scale of 0.235 arcseconds 
per pixel. The spectra cover the wavelength range from 
4790 A to 5850 A with 0.38 A per pixel and a median 
instrumental resolution of Oinst — 39.3 kms~^ (as mea- 
sured on the 5577 A night-sky line). The seeing varies 
from 1.2 to 4 arcseconds with a median value of 2.2 arc- 
seconds. Integration times vary from 1.800 s to 3.800 s 
and on-object exposures are typically split into two for 
cosmic rejection. For large galaxies where the DSS image 
of the galaxy exceeds the slit length, we obtain separate 
exposures of empty sky with an exposure time of 420 s 
at the end of the science observation. In order to avoid 
an azimuth move of the telescope, the sky exposures are 
typically obtained one hour in RA after the object but 
at similar DEC. 



Furthermore we observe a collection of kinematic tem- 
plate stars (G and K giants, see Tab. |3j metallicity: 

^ http://leda.uiiiv-lyoiil.fr 
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Fig. 1. — Distribution of Hubble types, central velocity dispersions and total magnitudes in the sample. 



[Fe/H] — -0.35 - 0.46) at the beginning and spectro- 
scopic standards throughout the duration of this cam- 
paign. The stars are wiggled and trailed along the slit 
such that a spectrum is recorded at each position where 
the star crosses the slit. This is used to map out the 
anamorphic distortion of the spectrograph. 



4. DATA REDUCTION 

We reduce the long-slit spectra following standard pro- 
cedures of bias subtraction, cosmic ray rejection and fiat 
fielding under MIDAS described in |Mehlert et al.| ( |2000" l 
with additional steps needed to correct tor spectral align- 
ment and anamorphism. We correct a two degree tilt 
between the spectra and the CCD rows by appropriate 
sub- pixel shifting of the CCD columns. Two bad columns 
at positions corresponding to A = 4850 A are corrected 
through interpolation. We perform the wavelength cali- 
bration on neon and cadmium arc frames with typically 
ten lines. Where the line signal is low we bin over a few 
rows along the spatial direction but never over more than 
five pixels corresponding to 1.2 arcseconds. After the 
original line identification we first fit a 4th-order poly- 
nomial to the line positions along the spatial direction 
in order to remove noise-induced row to row jitter, and 
then fit a 3rd-order polynomial along the spectral direc- 
tion. The remaining RMS scatter in line position is below 
one pixel. We then rebin the spectra in log-wavelength 
and correct for anamorphic distortion. The distortion of 
LRS is measured using stars which are trailed along the 
slit in order to generate several spectra or traces along 
the whole length of the slit. We centroid the traces by 
calculating the first moment of the photon count in a 
10 pixel wide window around the trace. We then first fit 
a 3rd-order polynomial to describe the trace position as 
function of wavelength and then a firther 3rd-order poly- 
nomial to the trace positions in each column to model 
the distortion. We find a distortion of up to ten pixel 
in the corners of the CCD with respect to the cen tre of 
the detector (see also Fig 2. in Sagha et al. 2010). We 
correct for the distortion by means of sub-pixel shifting. 
Counts of individual pixel are distributed into pixel of 
the target frame according the their overlapping surface 
area. We measure the distortion on several stellar spec- 
tra taken in similar manner at different nights. We find 
that the residual distortion at the edges of the chip — 
after correcting one stellar spectra with the distortion 



information of a different night — is never larger than 
1.5 arcseconds. This is below the typical FWHM of the 
PSF of our observations and, more importantly, well be- 
low the typical spatial bin sizes that we use at the ends 
of the slit. To correct for flexure of the instrument dur- 
ing the night we measure the wavelength position of the 
5577 A skyline at the slit ends and correct the wave- 
length calibration to zeroth order by adding a constant 
offset. The median absolute offset of all observations is 
17 kms"^. Where dedicated sky spectra are available, 
we collapse them along the spatial direction in order to 
obtain a single maximum signal-to-noise sky spectrum. 
This spectrum is then scaled according to the exposure 
time of the object and subtracted from the full frame. In 
cases where no sky frame is available, the sky signal is de- 
termined from the slit ends. One advantage of long-slit 
spectroscopy is that often the slit ends do contain suf- 
ficient non-object contaminated sky. However, the dif- 
ferential slit illumination is subject to change with time 
because the HET prime focus assembly moves across the 
telescope pupil during the duration of an observation. 
We test the effect of this differential illumination on 44 
blank sky spectra obtained over the course of this sur- 
vey. We use the slit ends to determine the sky signal in 
the same way as we do for the galaxy spectra. We then 
determine the differences between those and the sky sig- 
nal that we measure from the slit centre. We find that 
the residuals amount to no more than 5% of the sky sig- 
nal in all cases. We then derive kinematics using 5% 
higher and 5% lower sky values. The resulting errors 
are significantly smaller than the reported error bars in 
all cases. In the case of the major axis observation of 
NGC 3368, NGC 4569, and the minor axis observation 
of NGC 4569, the use of the dedicated sky frame results 
in an over-subtraction of the sky (i.e. clearly negative 
residuals) possibly because of stray light or an increased 
level of sky background at the time the sky frame was 
taken. In these cases we use the sky from the slit ends 
instead. In the cases of large galaxies such as NGC 3031, 
we test for object contamination by using different win- 
dow sizes at the slit ends for the sky extraction. We find 
the effect of object contamination to be negligible in all 
cases. Finally, all frames go through an extensive visual 
inspection. Artifacts like residuals of comic ray removal 
are corrected though interpolation of the neighbouring 
pixel. 
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TABLE 1 
Galaxy sample. 
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1) Galaxy name. 2) Hubble typo (RC3). 3) D 
Source for distance: l=p"ully (1994) 2= Tonry 

1991) 4= Saha 
Springob ct al 



et al 
([20091 



2001 


3= 


dc Vaucoulcurs ct al. 


2006 


5= 


tuUy ct al. (2009) 6= 



j|2009[ 



Blakcslcc c t al.[ 

9= |R ess eT^ l. (2009) 10= [Jensen ct al. ( 2003f~5) ' Total 
B-band magnitude (Hyperleda) . 6) Inclination (Hyper- 
Icda). Comments: a) No value in Hyperleda, from Peter 
Erwin (private communication). 

4.1. Derivation of the Kinematics and Template Library 
We derive stellar kinematics us ing the Fourier Correla- 



tion Quotien t (FCQ) method of Bender (19901; Bender 



et al.| ( |1994[ ). The log- wavelength cahbrated and sky 
subtracted spectra are spatially binned to reach S/N- 
values of at least 30 per pixel. An 8th-order polynomial 
is then fitted to the continuum and the first and last three 
channels in Fourier space are filtered out to remove low 
and high frequency variations in the continuum level. 

FCQ measures the full line of sight velocity distribu- 
tion (LOSVD). By means of deconvolving the autocorre- 
lation function, the FCQ method is more robust against 




4900 5000 5100 5200 5300 5400 
A [A] 



Fig. 2. — Continuum removed spectrum in the central radial bin 
of NGC 2841 (black) and the broadened template spectrum (red). 
Grey bars mark the positions of the [Olll] and [Nl] emission lines. 
Upper panel: Choosing the G8III, [Fe/H\ = 0.16 star HR3427 
results in a notable mismatch around the Mg triplet region. Best 
fit parameters are a = 235.1 ±3.0 kms"!, /13 = 0.037± 0.009, = 
0.041 ± 0.009, RMS = 0.034). Lower panel: Using HR2600 (K2III, 
[Fe/H] = —0.35) results in a much better match with cr = 241.2 it 
3.2 kms-i, /13 = 0.022 ± 0.009, hi = 0.048 ± 0.009, RMS = 0.024. 
While FCQ finds values for the LOSVD moments that fully agree 
within the errors, the residual spectrum will look very different for 
those two cases and render the detection of weak emission lines 
impossible in the case of HR3427. 

template mismatches than other Fourier or pixel-space 
based methods. Nevertheless, nebular emission lines can 
significantly affect the derived higher moments of the 
LOSVDs and therefore need to be taken into account 
for the derivation of the kinematics. This introduces the 
necessity of a very accurate model spectrum because oth- 
erwise residual mismatches between the observed galaxy 
spectrum and the broadened model spectrum will mimic 
emission signatures which are then incorrectly removed. 
We therefore form a pool of template spectra by com- 
bining actual observed stellar spectra (see Tab. Isl) with 
synthetic simple stellar population templates (SSP) from 
Vazdekis (1999). Those include varying metallicities and 
ages. We use a su b-sample of th e pubhshed spectra 
with Salpeter IMF (Salpeter [1955 ), and all combina- 
tions of ages of 1, 2, 5, 10, and 17.78 Gyr and metallic- 
ities of [Fe/H] = -1.68, -1.28, -0.68, -0.38, +0.00, +0.20. 
The published SEDs have a nominal resolution of 



1.8 A (FWHM) which corresponds to a* = 45 kms^^ 
and therefore slightly lower than the spectral resolution 
of dinst = 39.3 kms~^. 

We run FCQ with the collection of all velocity tem- 
plates. Then we choose the single best-fitting template 
based on the minimum RMS between the broadened tem- 
plate and the galaxy spectrum 



RMS = 



iG{X) - B{X) ® s,{x)y dX, 



(1) 



where Ai = 4817 A to A2 = 5443 A is the fitted wave- 
length range, G{x) is the galaxy spectrum, Si{x) is the 
i-th template spectrum, and B(x) is the broadening func- 
tion derived from FCQ. Note that this is different from 



Kinematics across bulge types 



5 



TABLE 2 
List of observations. 



Galaxy axis date seeing exp.-time angle sky 
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1 


800 


170 


no 


NGC 2880 


MJ 


2009- 


11- 


16 


3.4 


2 


200 


142 


no 


NGC 2880 


MN 


2009- 


12- 


18 


2.2 


2 


400 


52 


no 


NGC 2964 


MJ 


2010- 


02- 


18 


1.9 


2 


400 


96 


yes 


one; A 

NGC 2964 


MN 


2010- 


03- 


21 


2.2 


2 


400 


7 


yes 


NGC 3031 


MJ* 


2007- 


02- 


22 


2.2 


2 


700 


137 


yes 


NGC 3031 


MN'^ 


2005- 


12- 


28 


3.1 


1 


800 


67 


no 


NGC 3166 


MJ 


2008- 


02- 


06 


2.0 


2 


454 


85 


yes 


NGC 3166 


MN 


2008- 


12- 


25 


2.4 


2 


400 


175 


no 


NGC 3245 


MJ 


2008- 


02- 


06 


1.7 


2 


700 


174 


yes 


JNCjC o24o 


MN 


2008- 


12- 


25 


2.5 


2 


400 


84 


no 


NGC 3351 


MJ* 


2008- 


02- 


09 


1.5 


2 


550 


165 


yes 


M o o c; 1 

NGG 3351 


MN 


2008- 


12- 


27 


5.0 


2 


400 


75 


yes 


NGC 3368 


MJ* 


2007- 


02- 


26 


3.1 


2 


420 


153 


yes 


■NT /-I Q Q fl Q 


iVliN 


2008- 


12- 


no 


Z.D 


z 


4UU 


Q7 
O / 


yes 




MJ 


2009- 


12- 




1 7 
-1- . / 


2 




OO 


yes 




iViiN 


2010- 


02- 


1 Q 


1 S 
1 .o 


9 

z 


400 


i40 


yes 




A/T T 


2007- 


04- 


1 

lo 


1 .0 


9 

z 


700 




yes 






2009- 


01- 


uo 




9 
Z 


^9S 

ozo 


74 


yes 




A/T T 
iViJ 


2010- 


02- 


1 7 


1 9 
1 .Z 


9 
Z 


400 


H4 
04 


yes 


NGC 3dz / 


MJ 


2006- 


12- 


27 


2.3 


1 


800 


10 


no 






2007- 


02- 




9 9 


1 
± 


ftOO 


1 00 
iUU 


no 




A/T T 
iViJ 


2008- 


03- 


uo 


9 P, 
Z .D 


9 
Z 


700 


1 7S 


yes 




A/T T 
iViJ 


2007- 


04- 


ly 


1 .D 


9 
Z 


700 
/ UU 


1 OS 
iUo 


no 




Al Tfa 
iViJ 


2009- 


12- 


1 1 


9 1 
Z. i 


o 
z 


A nn 
4UU 


1 t; /I 
io4 


yes 


NGC 3945 


MN 


2010- 


04- 


12 


1.6 


4 


200 


64 


ves 


NGC 3953 


MJ* 


2008- 


02- 


06 


2.0 


2 


700 


32 


yes 


NGC 3992 


MJ* 


2008- 


12- 


28 


2.7 


2 


700 


66 


yes 


NGC 4030 


MJ 


2005- 


04- 


05 


2.3 


1 


800 


27 


no 


NGC 4203 


MJ 


2007- 


04- 


12 


1.3 


2 


520 


7 


yes 


NGC 4260 


MJ 


2008- 


12- 


29 


2.7 


2 


700 


62 


no 


NGC 4274 


MJ 


2007- 


04- 


19 


1.7 


2 


623 


99 


yes 


NGC 4314 


MJ 


2007- 


02- 


20 


2.3 


2 


700 


127 


no 


NGC 4371 


MJ* 


2006- 


12- 


27 


2.4 


1 


800 


85 


no 


NGC 4371 


MN 


2006- 


06- 


19 


1.8 


1 


800 


175 


no 


NGC 4379 


MJ 


2007- 


02- 


21 


2.7 


1 


800 


97 


yes 


NGC 4394 


MJ 


2007- 


05- 


11 


1.5 


2 


556 


123 


yes 


NGC 4448 


MJ* 


2007- 


04- 


16 


2.3 


2 


700 


85 


yes 



other algorithms such as t he Maximum Penalized Likcli- 

(j2000| or the 
ICappellari &: 



hood (MPL) technique of Gebhardt et al. 



Penalized Pixel-Fitting method (pP-X-i" ) of (Jappell 



Emsellem 
templates. 



) which fit a linear combination ot their 



subsequently fits G aussians with Her 
mite expansions (h^ and moments; Gerhard|1993 vai 
der Marel fc Franx|[T993l ) to the derived LUSVUs. In 
J^'ig.EJwe show examples of fits with two difi^erent broad- 
enedtemplates. While FCQ indeed finds very similar 
values for the moments of the LOSVD, the quality of the 
template match differs significantly in the two cases. 

In Fig. [3j we compare the impact of the usage of 
either just our observed templates or just the SSP li- 
brary. We compare only a subset of galaxies (NGC 2775, 
NGC 2880, NGC 3675, NGC 4030, NGC 4371 NGC 4379, 
and NGC 7457) for which we detect no significant emis- 



sion in order to assure that the derived moments are 
not affected by emission. The biases that we intro- 
duce by adding the SSP templates to our library are 
generally small (Act = —0.65 kms~^, A/13 = —0.011, 
A/14 = —0.012 and much smaller than our median errors 
on the respective moments. 

As our spectra often reach into the disk regions we deal 
with relatively low velocity dispersions. In a few cases the 
derived dispersions are of the order of the instrumental 
resolution. The matter gets complicated by the fact that 
the disk regions are also the regions of lowest surface 
brightness and therefore the regions with poorest S/N. 
It is important to understand how reliable the derived 
moments are under these circumstances. 

A caveat of the deconvolution in Fourier space is the 
amplification of high frequency noise. Fourier-based al- 
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TABLE 2 

— Continued 



Galaxy 


axis 


date 




(1) 


(2) 


(3) 




NGC4501 


MJ 


2010- 


04- 


06 


NGC4501 


MN 


2010- 


04- 


08 


NGC4536 


MJ 


2010- 


04- 


08 




MN 


2010- 


04- 


10 


JNCjC 45d9 


MJ 


2007- 


06- 


15 


iNLrL' 4oDy 


MiN 


2010- 


04- 


12 




A/T T 


2008- 


12- 


oc 

Zo 




A/T T 
iViJ 


2009- 


12- 


1 9 






2008- 


04- 


Ui 


NGC 4772 


MJ* 


2008- 


12- 


29 


NGC 4826 


Mj' 


2008- 


01- 


09 


NGC 4820 


MN'^ 


2009- 


06- 


25 


JNGG OU55 


MJ 


2008- 


03- 


05 


NGC 5055 


MN 


2009- 


06- 


26 


NGC 5248 


MJ' 


2007- 


04- 


18 


NGC 5566 


MJ 


2005- 


07- 


07 


NGC 5566 


MN 


2006- 


05- 


24 


NGC 7177 


MJ*- 


2007- 


08- 


11 


NGC 7177 


MN': 


2009- 


11- 


12 


NGC 7217 


MJ 


2007- 


08- 


12 


NGC 7217 


MN': 


2008- 


12- 


29 


NGC 7331 


MJ 


2007- 


08- 


11 


NGC 7743 


MJ*" 


2008- 


12- 


28 


NGC 7743 


MN': 


2009- 


10- 


17 



seeing exp.-time angle sky 



["] 




[°] 




(4) 


(5) 


(6) 


(7) 


2.0 


3,340 


140 


yes 


2.1 


2,505 


50 


yes 


3.6 


2,385 


120 


yes 


2.1 


2,500 


30 


yes 


O 1 
Z. i 


z, / UU 


1 A 

14 


yes 


a 


z,oou 


1 -LO 


yes 


2 4 


9 700 
Z, ( uu 


1 


no 


1.9 


2 400 


105 


yes 


9 9 


9 700 

z, / uu 


•^0 

ou 


yes 


9 


9 700 
Z, / UU 


140 


no 


2.0 


2,187 


96 


yes 


1.0 


9 1 99 
Z, IZZ 


9^ 
ZO 


yes 


4.0 


2,700 


103 


yes 


2.2 


2,400 


13 


yes 


1.6 


2,700 


109 


yes 


2.5 


1,800 


30 


no 


2.0 


1,750 


120 


no 


2.3 


2,700 


60 


yes 


1.9 


2,600 


173 


yes 


1.7 


2,700 


81 


yes 


1.4 


2,400 


178 


no 


1.7 


2,700 


171 


yes 


2.3 


2,400 


100 


no 


2.3 


2,400 


167 


no 



Notes: j^-j Q^jg^yry name. 2) MJ=major axis, MN=minor axis. 3) Date of 
observation. 4) Seeing FWHM. 5) Total exposure time. 6) Slit position 
angle east of north. 7) Dedicated sky frame was taken. 
Comments: a) No seeing information available, b) The position angle 
differs by more than 10° from the Hyperlcda published value, c) MN 
axis PA not orthogonal to MJ axis PA. 



TABLE 3 
Observed kinematic templates. 




MS 



-0.1 0.0 0.1 
[km/s] obs. tempi 




-200-100 100 200 
V [km/s] obs. tempi 



0.1 - 



0.0 - 



-0.1 - 




-0.1 0.0 0.1 
/14 [km/s] obs. tempi 




50 100 150 200 250 
a [km/s] obs. tempi 



Fig. 3. — Moments of the LOSVD that we obtain from the kine- 
matic extraction using the Vazdekis SSP library vs. the values that 
we obtain for the observed templates. The solid lines correspond 
to a one-to-one correlation, the dashed lines are actual fits to the 
data. In the case of v and a the fitted line is covered by the one- 
to-one line and not visible. Note: We only compare galaxies with 
no obvious sign of emission. 



Identifier 


type 


[Fe/H] 


date of obs. 


(1) 


(2) 


(3) 


(4) 


HR2600 


K2III 


-0.35 


04/03/2005 


HR3369 


G91II 


0.17 


04/02/2005 


HR3418 


K2III 


0.09 


04/03/2005 


HR3427 


KOIII 


0.16 


04/03/2005 


HR3428 


KOIII 


0.24 


04/03/2005 


HR3905 


K2IIIb 


0.46 


04/02/2005 


HR6018 


KOIII-IV 


0.01 


04/02/2005 


HR6159 


K7III 


-0.13 


04/02/2005 


HR6770 


G81II 


-0.05 


04/03/2005 


HR6817 


KlIII 


-0.06 


04/02/2005 


HR7576 


K3I11 


0.42 


04/03/2005 



Notes: Identifier. 2) Stellar classification". 

3) Metallicity*^ . 4) D ate of observation. 
Comments: a) From IWorthey et al.| | [T994l l. 

gorithms therefore filter tlie signal before the actual de- 
convolution step. FCQ uses the o ptimal Wiener filter 
(iBrault & White|1971l|Simkin|1974|). The basic idea is to 



decompose the Fourier transform of the input data into 
a Gaussian contribution — the data part - and an expo- 
nential function — the noise part. The optimal Wiener 
filter then weighs the various signal channels according 
to their relative contribut ion to the da ta part of the in- 
put signal (for details see Bender|["l990 1. While a purely 
Gaussian LOSVDs ought to be well modelled by a Gaus- 
sian in Fourier space, the Gauss-Hermite moments cause 
higher frequency shoulders, that are easily swallowed by 
the noise. Adjusting the filter width may rec over charac- 
teristics of the LOSVD (iBender et al. 11 19941) at the cost 
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TABLE 4 
Parameters for the linear bias 
corrections in velocity dispersion 



TABLE 5 

Parameters for the linear bias corrections 
IN /i3 and /l4 







S/N per 


pixel 






22.5 


40.0 


37.5 


75.0 




1.06 


1.06 


1.06 


1.07 




-11.04 


-10.89 


10.80 


-10.70 



of increased statistical uncertainty. 

Here we choose not to broaden the Wiener Filter as this 
yields better stability against statistical deviations. But 
this causes biases especially at low velocity dispersions. 
To correct for these biases we carry out extensive Monte 
Carlo simulations on a regular parameter grid of velocity 
dispersion, h^, h^, S/N and template. We generate arti- 
ficially broadened spectra at each grid point with 30 dif- 
ferent noise realisations according to the input signal to 
noise. We find that the necessary corrections to a, and 
/i4 are well behaved and linear functions between input 
and retrieved values and independent of input template if 
the velocity dispersion is larger than 75 kms~^, the signal 
to noise is larger than 30 per pixel and a stellar template 
is used. The SSP templates cause non-linear behaviour 
at small velocity dispersions. While we still use the SSP 
templates to generate broadened model spectra during 
the emission line fitting, the reported kinematic values 
are exclusively based on stellar templates, and corrected 
for biases using 

Tab. [2] and [5l list the corresponding parameters that we 
obtain from the simulations. For velocity dispersions be- 
low 75 kms^^ and S/N < 30 per pixel, we do not report 
values for and h^. Further we report values of v and 
a only for S/N > 20. 
We estimate statistical errors in the derived moments 



throu g h Mon te Carlo simulations as described in Mehlert 
et al.| (|2000|). Once the optimum LOSVD is derived by 



FCnTwe generate the synthetic spectra using the fitted 
V, a, /13 and /14-parameters and the best fitting stellar 
template. In a similar manner as for the bias correction, 
100 different realisations of artificial noise are then added 
to the spectra to reach the same signal to noise values 
as in the original spectra. We then use FCQ again to 
derive the kinematics on those spectra. The reported 
errors correspond to the statistical one-sigma deviations 
from the mean. 



4.2. Emission line subtraction and gas kinematics 

A significant fraction of the objects in our 
sample show emission in Hp (4861.32 A), 
[Oiii] (4958.83 A & 5006.77 A) and [Ni] 
(5197.90 A & 5200.39 A). The nitrogen emission 
line lies on the red flank of the Mg triplet feature - 
the most important kinematic feature in our spectral 
range. While typically weak, the nitrogen emission 



cr[kn 



S/N per pixel 
30.3 37.5 75.0 





75.0 


1.3084 


1.2947 


1.2734 




100.0 


1. 1281 


1 . 1 142 


1 .0874 




150.0 


1.0182 


1.0104 


1.0000 




200.0 


1.0103 


1.0059 


0.9988 




250.0 


1.0037 


1.0008 


0.9945 




75.0 


0.0002 


0.0001 


0.0001 




100.0 


0.0003 


0.0003 


0.0002 




150.0 


0.0003 


0.0003 


0.0003 




200.0 


0.0009 


0.0008 


0.0006 




250.0 


0.0010 


0.0010 


0.0009 




75.0 


1.8521 


1.8088 


1.7277 




100.0 


1.4994 


1.4500 


1.3555 




150.0 


1.0857 


1.0655 


1.0307 




200.0 


1.0407 


1.0275 


1.0045 




250.0 


1.1503 


1.1371 


1.1089 




75.0 


0.0864 


0.0814 


0.0738 




100.0 


0.0280 


0.0240 


0.0160 




150.0 


0.0165 


0.0148 


0.0127 




200.0 


0.0105 


0.0104 


0.0104 




250.0 


-0.0027 


-0.0015 


0.0006 



often significantly affects the derivation of ft,3-moments. 
ft,3-moments measure the asymmetric deviation from a 
Gaussian and are expected to behave antisymmetrically 
with respect to the galaxy centre in the case of axisym- 
metric systems. Deviations from this antisymmetry may 
hint at contamination by nitrogen emission lines. We 
therefore decided to remove nebular emission following a 
simila r procedure as the GANDALF routine ( Sarzi et al. 
2006 1: We perform a first fit to the galaxy spectrum 



over a larger spectral range reaching from 4820 A to 
5440 A using the FCQ algorithm. We then subtract the 
best fitting broadened stellar spectrum from the galaxy 
spectrum and fit Gaussian functions — using a standard 
least squares algorithm — to the residual emission. The 
algorithm first searches for emission in a 500 kms~^ 
window around the brighter oxygen line at 5007 A (the 
oxygen doublet is well resolved at our instrumental 
resolution) red-shifted to the systemic velocity. It fits 
for the three parameters of amplitude, central velocity, 
and the velocity dispersion. It then goes on to the other 
and generally weaker emission lines and performs a fit to 
their amplitude while assuming the same central velocity 
and velocity dispersion as the oxygen line. In principle 
the ratio of the two oxygen emission lines is given by 
atomic physics and is a constant of value 0.33. Rather 
than fixing these values during the fit we also fitted the 
lower amplitude line as this provided another handle 
on the reliability of our method. We then subtract 
the best fitting emission lines from the original galaxy 
spectrum and repeat the FCQ multiple-template fit. 
The best fitting broadened template is again subtracted 
from the input spectrum and the gas emission fit is 
repeated on the improved difference spectrum. An 
example for a spectrum that shows signs of nebular 
emission is shown in the upper panel of Fig. [4] We 
plot the residuals between the recorded spectrum and 
the broadened model spectrum after the removal of the 
emission in the lower panel. This iterative approach 
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4900 5000 



5100 5200 
A [A] 



5300 5400 



HST F438W, but it is misaligned with the outer disk also 
seen only in vicinity of the nucleus. 

Yet, a few objects remain for which we do not 
feel confident to assign a classification based on their 
HST morphology: NGC2460, NGC3953, NGC4826 and 
NGC 7217. We treat them as unclassified throughout 
this work. As mentioned in fj2j we also do not classify 
the bulges of the galaxies NGCJ 3593 and NGC 7331 due 
to their high inclination. In Appendix IB] we give a de- 
tailed explanation for the bulge classincation for each 
individual object. 

4.4. Photometry 

We use the results from decompositions of surface 
brightness profiles to investigate possible correlations be- 
tween photometric parameters and kinematic structure. 
Also the bulge-disk decompositions serve to identify the 
actual bulge region of a particular galaxy. 



Fig. 4. — Example of the nebular emission line subtraction. Here 
the gas emission was removed from one of the central spectra of 
NGC 3368. Upper panel: After subtraction of the best fitting 
broadened stellar template (red) the algorithm finds a significant 
emission signal in the [Oni] lines and the nitrogen doublet (green). 
The vertical bars mark the search range for emission. Lower panel: 
Residuals between observed and model spectrum after removal of 
the emission. 

was found to converge very quickly. A fourth FCQ 
fit typically yields no significant change in the derived 
kinematics any more. Our reported stellar kinematics 
went through three subsequent iterations of template 
fitting with two interleaved gas emission removal steps. 
The gas re moval routine give s access to the study of 



line ratios (Sarzi et al. 2006) and is also a necessary 
step before the calculation ot absorption line indices and 
subsequent stellar population analysis. Both of which 
will be subject to a forthcoming publication. 

4.3. Identification of pseudobulges 

If no classi fica tion is available already from [Fisher fc] 
roryl ( |2008l ) or [Fisher fc P rory' ('2010[), we follow the 
same procedure tor the identihcation of pseudobulges. 
We define a bulge photometrically as the excess light 
over the inwards extrapolation of the outer disk expo- 
nential luminosity profile. Th e bulge-to-disk decompo- 
sitions that we adopt in §4.4[ allow us to determine the 
bulge region of an object. Here, we classify bulges us- 
ing close-to y-band HST images (i^547M, F555W & 
F606W). While these bands are subject to dust obscu- 
ration they are also sensitive to an enhancement in star 
formati on rate, a fea ture commonly observed in pseudo- 
bulges (Fisher 2006). We visually inspected the HST 
images to see whether the bulge regions contain disk-like 
structures such a nuclear spirals, nuclear bars and/or nu- 
clear rings. If such structure is present we call the pho- 
tometric bulge a pseudo-bulge. If there is no structure 
(the bulge resembles an elliptical galaxy with a smooth 
light distribution), we call this bulge a classical bulge. 

Weak central dust s pirals that also occur frequently in 
elliptical galaxies (e.g. Storchi-Bergmann et al.[[2007[) — 
often distinguishable from the outer disk because they 
are inclined differently — are no reason for us to call a 
bulge a pseudobulge. NGC 2841 poses an example for 
this situation. A dust spiral can easily be identified in 



We de rive surface photo metry follow ing the prescrip 



tions m 



Fisher fc Drory[ p008j and [Fisher fc Drory 



([2010|. I'br each galaxy, we combine multiple data 
sources mostly in the infrared (but sometimes in optical 
bands) to obtain a final ID composite surface bright- 
ness profile. The different data have been calibrated 
against the H-band using 2MASS magnitudes. High- 
resolution HST imaging is used in the galaxy center, 
while wide-field images sample the outer disk. The re- 
sulting composite profile of each galaxy is used to de- 
rive the bulge-to-disk photometric decomposition. Our 
method is well tested and has been used in several pub 



licati ons (Fisher & Drory 2008 



2009 1 . Our principal sourc e of data is 2MASS H-band 



2010 Kormendy et al. 



maps ( Skrutskie et al. 2006 1 . When available we use dat a 



from 2MASS Large Galaxy Atlas (Jarrett et al. 2003) 



We prefer iJ-band over because the 2MASS if -band 
data is more sensitive than the 2MASS Ks data. For all 
galaxies we use the NASA/IPAC Extragalactic Database 
to search for ancillary ground based iJ-band data. We 
also include Spitzer 3.6 /im data. Finally, when avail- 
able we also include high resolution F160W images from 
HST/NICMOS. In a few cases, where galaxies were lack- 
ing archival NICMOS data, we use I or R band data 
instead. The high resolution data can be crucial t o accu- 
rately con straining the bulge-disk decomposition. |Fisher| 
fc Drory] ([2010) investigates the uncertainty introduced 
from mixing filters in this way, it is typically smaller 
than 0.1 mags, and there fore small compared to the un- 
certainty in the fit. Also, Fisher & Drory (20081 derives 
very similar Sersic indices with V-band data as with H- 
band profiles. 

We fit ellipses to all i mages. Isophotal fitting is car - 
ried out using the code of Bender fc Moellenhoff (1987). 
See [Fisher fc Droryl ( [2008| ) for a brief summary of the 
procedure. I'he code returns a 2-D surface brightness 
profile (including for each ellipse center, major fc minor 
axis size, position angle, and mean surface brightness). 
We then combine all profiles into a composite surface 
brightness profile. The power of this method is two-fold. 
First, combination of surface brightness profiles allows us 
to robustly identify systematic errors from point-spread- 
functions and sky subtraction. Secondly the resulting 
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TABLE 6 
Data Sources 



No. 


IiistrumGiit 


Filter 


Scale 


Field Of View 








(arcsec/pixcl) 


(arcscc) 


1 


HST Nicmos 1 


F160W 


0.043 


11 X 11 


2 


HST Nicmos 2 


F160W 


0.075 


19.2 X 19.2 


3 


HST Nicmos 3 


F160W 


0.20 


51.2 X 51.2 


4 


2MASS 


H 


1.0 


512 X 1024 


5 


2MASS 


H 


1.0 


variable size, 










mosaic 


6 


Spitzcr IRAC 


3.6 /im 


0.60 


variable size, 










scan 


7 


Spitzcr IRAC (SINGS) 


3.6 /im 


0.75 


variable size. 










scan 


8 


HST ACS/WFC 


F814W 


0.049 


202 X 202 


9 


HST WFPC2 


F547M 


0.10 


80 X 80 


10 


HST WFPC2 


F814W 


0.10 


80 X 80 


11 


Perkins 1.8 m OSIRIS 


H 


1.5 


412 X 412 


12 


Lick 3 m pNIC 


K 


0.24 


15.4 X 15.4 


13 


CTIO 1.5 m OSIRIS 


H 


1.1 


312 X 312 


14 


William Hcrschcl Telescope 


K 


0.24 


252 X 252 




INGRID 








15 


Calar Alto Observatory 2.2 m 


K 


0.66 


172 X 172 




MAGIC Nicmos 3 








16 


Mauna Kea 0.61 m Nicmos256 


K 


2.1 


644 X 568 


17 


UKIRT 3.8 m IRCAM II 


H 


1.7 


198 X 72 



Reference 



MAST archive^ 
MAST archive^ 
MAST archive^ 
IPAC archive 



Jarrett et al. 2003 



IPAC archive 



SINGS'" 



MAST archive^ 
MAST archive^ 

MAST archive- ^ 

[Eskridse et al. V2OO2I 
Hausctier (1995) 
Eskridge et al.^l,2002[ 
Knapeii et al.lJzDD^lr 



MoUenhott &; Heidt 2001 



l^lly et al.|||19ybf 



g Hz varTaer Kruit| | |1994[ 



1) http; //archive . stsci . edu/ 2) http : //irsa. ipac . caltech . edu/ 3) http : //data . spitzer . caltech . edu/popular/sings/ 



composite profile has an extremely high dynamic range 
in radius, which is necessary to accurately constrain the 
b ulge-disk decomp osition (see dis cussion in the Appendi x 
of|Fisher & Drory||2008| and also |Kormendy et alj2009|). 



The zero points ot our profiles are matched against the 
2MASS data. 

We determine bulge and disk parameters by fitting 
each surface brightness profile with a one-dimensional 
Sersic function plus an exponential outer disk, 



/(r) = /oexp \-{r/royH + exp h(r//i)] 



(2) 



where r represents the distance along the major axis, Iq 
and To are the central surface brightness and scale length 
of the bulge, Id and h represent the central surface bright- 
ness and scale length of the ou ter disk, and n represents 
the bulge Sersic index (Sersic 1968| ). The half-light ra- 
dius, Te, of the bulge is obtained by converting tq. 



(3) 



where the value of 6„ is a proportionality cons tant de- 
fined such that r(2n) = 27(2n,6„) (Ciottil |199l"| ). F and 
7 are the complete and incomplete gamma functions, re- 
spectively. We use the approximation 6„ « 2.17n — 0.355 
(Caon et al. 1993). We restrict our range in possible 
Sersic indices to n > 0.33 to ensure that the approxima- 
tion is accurate. Bulge and disk magnitudes are adjusted 
to account for the shape of the bulge using the ellipticity 
profile from the isophote fitting. 

Intermediate type galaxies are known to contain many 
components that are not well described by the decompo- 
sition into a Sersic bulge and exponential disk (e.g. bars, 
rings, nuclear star clusters). Similar to the outer disk, in 
the bulge we exclude significant non-Scrsic components 
such as nuclear bars and nuclear rings. The inclusion 
of these features can have unpredictable effects on the 
Sersic index, depending on the relative size of the fea- 



will affect the decomposition. Essentially, this has the 
effect of decreasing the robustness of the fit, which will 
be refiected in the error bars. The parameters from the 
decompositions are presented in Tab. [7] , the image data 
sources are described in Tab. |6l 

The surface brightness profile of NGC 5566 does not 
follow a typical bulge/disk profile. We publish the ob- 
tained kinematics here but exclude this object from all 
further analysis. 



4.5. Bulge Radius 

Here we are particularly interested in the kinematic 
properties of the bulge regions of our observed galaxies. 
Of course the derived LOSVDs will always be the light 
weighted average of all components (bulge, disk, bar) 
along a particular line of sight through the galaxy. But 
the photometric bulge to disk decompositions allow us 
to determine within which radius the bulge should dom- 
inate. We define the bulge radius along the major axis 
as the radius where the light contribution of the photo- 
metric bulge component exceeds the light contribution of 
the disk component by 25%: 



bulge 



exp 



1.25 ■ I^''" exp 



(4) 



where /^"'^"^ is the central surface brightness of the bulge 
component, J^***^ is the central surface brightness of the 
disk component, tq and h are the scale lengths of the 
bulge and disk components and n is the Sersic index (see 
previous Section for the relation between tq and Tg). One 
might argue that the bulge effective radius is a more 
natural choice as Tf, of course is dependent on the disk 
parameters. But we find that in a number of galaxies 
actually lies in a region that is dominated by disk light. 
The choice of 25% is a compromise between the desire 
to be reasonably dominated by the bulge component on 



ture, and wh at type it is. The appendix of Fisher & 



Drory (2008) discusses how masking data in the bulge 
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TABLE 7 

Bulge to disk decomposition parameters". 



Galaxy 

(1) 


Bulge 
morph. 
(2) 






n 

(3) 






5 arcsec 
(4) 


-2 


arcsec 
(5) 


mag 
(6) 






maj 


disk 

Mo 

5 arcsec 
(7) 


-2 


h 

arcsec 
(8) 






771''' 

ma 
(9) 






Data sources 
(10) 


NGC 1023 


c 


2 


52 


± 





81 


15 


76 


± 


1 


11 


12.35 ± 4.39 


7 


15 


± 


1 


11 


16 


62 


± 





32 


62.20 ± 


6 


56 


5 


66 


± 





35 


5 


10,15 


NGC2460 


P 


3 


49 


± 





32 


18 


02 


± 





46 


12.69 ± 4.40 


9 


19 


± 





46 


16 


40 


± 





28 


11.42 ± 





98 


9 


11 


± 





31 


2 


4 


NGC 2681 


P 


3 


82 


± 





31 


14 


58 


± 





63 


3.74 ± 3.00 


8 


35 


± 





63 


17 


32 


± 





33 


23.78 ± 


2 


76 


8 


44 


± 





37 


3 


4,6 


NGC 2775 


c 


3 


23 


± 





93 


17 


28 


± 


1 


02 


15.86 ± 5.50 


8 


00 


± 


1 


02 


17 


20 


± 





61 


41.28 ± 


7 


91 


7 


12 


± 





67 


4 


6,10,11,14,15 


NGC 2841 


c 


3 


22 


± 





58 


16 


55 


± 





68 


15.46 ± 8.98 


7 


33 


± 





68 


16 


49 


± 





17 


60.51 ± 


2 


82 


5 


58 


± 





18 


2 


5,7 


NGC 2859 


c 


2 


34 


± 





65 


16 


23 


± 


1 


08 


8.3(< 21.4'') 


8 


52 


± 


1 


08 


19 


11 


± 





47 


55.17 ± 


8 


93 


8 


41 


± 





52 


5 


6,8 


NGC 2880 


c 


3 


41 


± 





48 


17 


49 


± 





59 


11.7(< 24.2'') 


8 


85 


± 





59 


18 


29 


± 





40 


25.47 ± 


1 


88 


9 


26 


± 





41 


4 


10 


NGC 2964 


p 


1 


01 


± 





34 


15 


43 


± 





51 


2.04 ± 0.50 


11 


18 


± 





51 


16 


40 


± 





18 


16.07 ± 





65 


8 


37 


± 





19 


2 


4,6 


NGC 3031 


c 


4 


09 


± 





48 


17 


14 


± 





62 


70.70 ± 54.91 


4 


49 


± 





62 


16 


59 


± 





25 


132.80 ± 


8 


05 


3 


98 


± 





26 


2 


5,7 


NGC 3166 


p 


1 


24 


± 





30 


14 


37 


± 





39 


4.36 ± 1.12 


8 


37 


± 





39 


15 


86 


± 





44 


15.47 ± 


2 


66 


7 


92 


± 





50 


5 


6,9,11 


NGC 3245 


c 


2 


75 


± 





56 


15 


16 


± 





80 


4.51 ± 1.70 


8 


69 


± 





80 


16 


44 


± 





28 


21.60 ± 


1 


49 


7 


77 


± 





30 


2 


4,6 


NGC 3351 


p 


1 


38 


± 





74 


15 


99 


± 





60 


8.08 ± 3.10 


8 


59 


± 





60 


17 


01 


± 





44 


49.48 ± 


5 


13 


6 


54 


± 





47 


2 


5,7 


NGC 3368 


p 


2 


46 


± 





77 


15 


97 


± 





75 


13.08 ± 7.29 


7 


25 


± 





75 


16 


57 


± 


1 


95 


35.26 ± 29 


61 


6 


84 


± 


2 


31 


2 


4,7 


NGC 3384 


p 


1 


58 


± 





22 


14 


39 


± 





40 


5.35 ± 1.16 


7 


83 


± 





40 


16 


96 


± 





15 


44.42 ± 


3 


13 


6 


73 


± 





18 


2 


5,6 


NGC 3521 


c 


3 


66 


± 





77 


15 


48 


± 


1 


50 


8.55 ± 6.36 


7 


48 


± 


1 


50 


15 


99 


± 





31 


49.53 ± 


3 


87 


5 


52 


± 





33 


2 


5,7 


NGC 3593 


7 


1 


22 


± 





21 


16 


18 


± 





27 


14.55 ± 2.16 


7 


57 


± 





27 


17 


61 


± 





30 


52.52 ± 


4 


99 


7 


02 


± 





33 


3 


5,6,11 


NGC 3627 


P 


1 


50 


± 





58 


14 


53 


± 





59 


4.77 ± 1.52 


8 


24 


± 





59 


16 


72 


± 





15 


65.99 ± 





15 


5 


63 


± 





16 


3 


5,7 


NGC 3675 


P 


1 


57 


± 


1 


12 


16 


35 


± 


1 


93 


9.0(< 23.4'') 


8 


66 


± 


1 


93 


16 


22 


± 





31 


36.85 ± 


3 


42 


6 


40 


± 





34 


2 


5,6,11 


NGC 3898 


c 


3 


22 


± 





86 


15 


93 


± 


1 


27 


7.63 ± 2.88 


8 


24 


± 


1 


27 


17 


25 


± 





64 


29.00 ± 


4 


56 


7 


94 


± 





68 


2 


5,6 


NGC 3945 


p 


1 


79 


± 





48 


16 


00 


± 





76 


9.82 ± 3.50 


8 


05 


± 





76 


19 


02 


± 





44 


83.17 ± 26 


76 


7 


43 


± 





64 


5 


6,10 


NGC 3953 


7 


2 


43 


± 





68 


17 


15 


± 





96 


12.74 ± 8.03 


8 


49 


± 





96 


17 


44 


± 





15 


66.14 ± 


4 


73 


6 


34 


± 





19 


2 


3,5,7 


NGC 3992 


C 


3 


18 


± 


1 


18 


17 


44 


± 


1 


44 


12.23 ± 3.74 


8 


73 


± 


1 


44 


17 


65 


± 





50 


77.54 ± 19 


36 


6 


20 


± 





62 


5 


9,16 


NGC 4030 


P 


1 


98 


± 


1 


30 


16 


50 


± 


1 


52 


5.18 ± 2.14 


9 


89 


± 


1 


52 


15 


60 


± 





36 


15.82 ± 


1 


95 


7 


61 


± 





41 


2 


4,6,13 


NGC 4203 


c 


2 


45 


± 





83 


15 


72 


± 


1 


51 


6.78(< 15.9'') 


8 


43 


± 


1 


51 


17 


28 


± 





37 


30.31 ± 


3 


53 


7 


88 


± 





41 


4 


6,10 


NGC 4260 


c 


3 


68 


± 





42 


19 


08 


± 





48 


21.49 ± 11.93 


9 


08 


± 





48 


17 


02 


± 





19 


21.69 ± 


1 


59 


8 


34 


± 





22 


2 


4,6 


NGC 4274 


p 


1 


52 


± 





24 


15 


49 


± 





28 


5.92 ± 1.10 


8 


73 


± 





28 


17 


03 


± 





16 


46.53 ± 


2 


92 


6 


70 


± 





19 


5 


6,8 


NGC 4314 


p 


2 


72 


± 





96 


17 


12 


± 


1 


09 


10.53 ± 6.96 


8 


82 


± 


1 


09 


16 


70 


± 





36 


35.08 ± 


2 


65 


6 


98 


± 





38 


2 


4,6,14 


NGC 4371 


p 


2 


21 


± 


1 


00 


16 


66 


± 


1 


35 


10.97 ± 6.00 


8 


37 


± 


1 


35 


18 


10 


± 





98 


44.59 ± 13 


45 


7 


86 


± 


1 


07 


1 


4,6 


NGC 4379 


c 


2 


39 


± 





55 


16 


72 


± 





69 


6.40 ± 1.97 


9 


56 


± 





69 


16 


91 


± 





37 


13.52 ± 


1 


48 


9 


26 


± 





40 


1 


4,8 


NGC 4394 


p 


1 


58 


± 





67 


16 


28 


± 





87 


6.10 ± 1.81 


9 


43 


± 





87 


18 


40 


± 





25 


57.75 ± 


6 


38 


7 


60 


± 





30 


4 


6,11,12 


NGC 4448 


p 


1 


19 


± 





25 


16 


43 


± 





31 


6.70 ± 1.01 


9 


52 


± 





31 


16 


85 


± 





17 


28.42 ± 





97 


7 


59 


± 





18 


5 


6,10,11 


NGC 4501 


p 


1 


25 


± 


1 


06 


15 


43 


± 


1 


34 


5.31 ± 1.53 


8 


99 


± 


1 


34 


15 


81 


± 





46 


39.26 ± 


3 


89 


5 


85 


± 





49 


5 


6,12 


NGC 4536 


p 


1 


47 


± 





35 


14 


77 


± 





62 


3.98 ± 1.18 


8 


88 


± 





62 


17 


18 


± 





15 


32.32 ± 


1 


82 


7 


64 


± 





17 


2 


5 


NGC 4569 


p 


2 


34 


± 





97 


15 


13 


± 


1 


59 


4.80 ± 2.33 


8 


60 


± 


1 


59 


16 


84 


± 





32 


61.32 ± 


6 


42 


5 


91 


± 





36 


2 


3,5,7 


NGC 4698 


c 


2 


51 


± 





53 


15 


66 


± 





66 


5.11 ± 1.41 


8 


97 


± 





66 


17 


32 


± 





34 


34.93 ± 


2 


76 


7 


61 


± 





36 


2 


4,6,11 


NGC 4736 


p 


1 


23 


± 





30 


13 


80 


± 





41 


7.76 ± 1.65 


6 


56 


± 





41 


14 


93 


± 





35 


26.82 ± 


3 


85 


5 


80 


± 





41 


3 


5,7 


NGC 4772 


c 


3 


03 


± 





88 


17 


73 


± 


1 


20 


13.49 ± 8.24 


8 


84 


± 


1 


20 


18 


87 


± 





55 


71.96 ± 16 


04 


7 


59 


± 





64 


4 


11 


NGC 4826 


7 


3 


93 


± 





88 


16 


93 


± 


1 


06 


28.93 ± 10.26 


6 


25 


± 


1 


06 


16 


50 


± 





27 


67.95 ± 


6 


86 


5 


34 


± 





31 


2 


5,7 


NGC 5055 


P 


1 


71 


± 


1 


03 


17 


35 


± 


1 


36 


27.88 ± 16.08 


7 


16 


± 


1 


36 


16 


57 


± 





38 


68.70 ± 


5 


53 


5 


40 


± 





40 


3 


5,7 


NGC 5248 


P 


1 


29 


± 





45 


16 


47 


± 





67 


0.7(< 4.9'') 


8 


74 


± 





67 


17 


43 


± 





28 


43.99 ± 


3 


72 


7 


21 


± 





30 


2 


5,6 


NGC 7177 


P 


2 


03 


± 





52 


16 


35 


± 





59 


7.02 ± 3.17 


9 


08 


± 





59 


16 


42 


± 





31 


15.45 ± 





73 


8 


49 


± 





32 


2 


4,6,17 


NGC 7217 


7 


3 


20 


± 


1 


03 


17 


04 


± 


1 


21 


13.41 ± 9.35 


8 


13 


± 


1 


21 


16 


50 


± 





64 


28.54 ± 


6 


60 


7 


23 


± 





72 


2 


4,6,11 


NGC 7331 


7 


2 


85 


± 


1 


02 


16 


20 


± 


1 


30 


16.63 ± 9.90 


6 


88 


± 


1 


30 


16 


97 


± 





47 


61.61 ± 


9 


04 


6 


03 


± 





52 


3 


5,7 


NGC 7743 


P 


3 


66 


± 





52 


15 


15 


± 





89 


2.37 ± 1.04 


9 


94 


± 





89 


17 


14 


± 





22 


22.27 ± 


1 


96 


8 


41 


± 





25 


2 


4,6 



° 1) Galaxy name. 2) Bulge classification: c — classical, p — pscudobulge, ? — not classified. 3) Bulge Scrsic index. 4) Bulge surface briglitness 
at Te- 5) Bulge efl^ectivc radius along the major axis. 6) Bulge apparent magnitude. 7) Disk central surface brightness. 8) Disk scale length. 9) Disk 
apparent magnitude. 10) Image data sources, see Tab .[6] 

Co7n77ie7its: a) These photometric decompositions arc Dased on infrared and optical data, but calibrated against the H-band. b) The error on the 
effective radius is comparable to or larger than the value itself. We list the value that is preferred by the fit and the upper limit. 



the one hand and still wanting to maintain a sufficient 
number of resolution elements within the bulge radius on 
the other. The values for the bulge radius are listed in 
Tab. [9j In Fig. [16] we indicate the location of the bulge 
radius though a oashed vertical line. 

5. RESULTS 

5.1. Kinematic profiles and comparison with literature 

Tab. |8] gives an example of the format of the measured 
stellar kinematic moments as function of the distance 
from the center of the galaxy. The full listing is available 
electronically!**] In the Appendix |X1 we plot the kine- 
matic profiles. When available, we also plot data from 
the literature for comparison. Integral Field Spectro- 
scopic data from SAURON is available for some of the 
galaxies in our sample. In those cases we create pseudo 

** http://cds.u-strasbg.fr/ 



long-slit data through interpolation of the SAURON v, 
cr, h^, and maps along a slit aperture with a posi- 
tion angle corresponding to our observation. In general 
the agreement of our data with the published values is 
acceptable. 

In a few cases such as NGC 4203 a difference between 
the previously published data and ours are explained by 
the difference in the observed position angle. 

B ertola et al.| (|1995 ) find somewhat larger velocity dis- 
persions tor JN GC 4379 than we do. Formally their instru- 
mental dispersion should allow to resolve the 80 kms^^ 
118 kms~^ that we find for the dispersion in the bulge. 



Dumas et al. ( 2007 ) finds larger velocity dispersions in 
the cases of iNGCJS^l and NGC 5248 than we do. The 
dispersion of those objects is probably too low to be re- 
s olved by their instrumental d ispersion of « 110. kms~^. 

(|2001j) find systematically lower 
/ 2841 than we do. However, 



Vega Beltran et al 
velocity dispersions tor 



Kinematics across bulge types 
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we also plot data from Heraudeau & Simien| (|1998|) which 
are in excellent agreement with ours 



The SAU RON data for NGC4698 ( |Falc6n-Barroso 
et al. 2006) suggest a somewhat larger velocity disper- 
sion over our whole observed range than we find. They 
also find negative /13 moments on the east side. The 
dispersion of this galaxy 140 kms""'^) should be well 
resolved by SAURON and such the difference remains 
somewhat mysterious but small. 



5.2. Signatures of bars in velocity profiles 

In our sample, 29 out of 45 of the galaxies are classi- 
fied as barred or as hosting an oval. Bars and ovals will 
affect the observed kinematics and their presence should 
be reflected in the moments of the observed LOSVD. 



Bureau &: Athanassoula ( 2005 1 use iV-body simulations 
to derive diagnostics for the presence of bars in edge- 
on disks. They find that double-hump rotation curves, 
plateaus and shoulders in velocity dispersion, and cor- 
relation of /13 moments with velocity in contrast to the 
usually-seen anti-correlation are indicators for the pres- 
ence of a bar. The double-hump describes a rotation 
curve that first rises quickly with radius, reaches a local 
maximum then drops slightly and starts rising again to- 
wards larger radii. We do see similar features in a number 
of our galaxies even though they are not observed edge- 
on (e.g. see the rotation curves for NGC 2841, NGC3351, 
and NGC 3384 in Fig. 16). The signature is not always 



strong enough to form an actual local minimum after 
the fast inner rise. Instead, in some cases we observe 
shelves: the rise in velocity stagnates for a certain radial 
range but becomes larger again before fin ally flattening 
out (e.g. NGC 1023 and NGC 3627 in Fig. [l6|. 

Out of 29 barred galaxies (including 6 ovals), 20 do 
show such features. However our data do not extend very 
far into the disk region in many of the objects in our sam- 
ple; also visibility may be inhibited by the coarse spatial 
binning of some of our spectra. Further, this diagnos- 
tic tool was developed for edge-on systems, so it is likely 
that we miss bar signatures in the velocity profiles. How- 
ever, 9 of the 16 non-barred galaxies show either shelves 
or double humps which may be an indication that those 
systems actually do host a bar that is not readily seen 
photometrically. 

5.3. Central velocity dispersions 

We calculate the central velocity dispersion of the 
galaxies in our sample by averaging the major axis dis- 
persion within a tenth of the effective bulge radius r^ that 
we obtain from the photometric decomposition. The val- 
ues for the central dispersions are given in Tab. |9| The 
quoted errors correspond to the formal errors of the de- 
rived mean within re/10. 

In Fig. [5] we show corresponding histograms of the cen- 
tral dispersions. In the left panel we discriminate bulge 
types based on their morphology, in the right panel we 
discriminate by Sersic index. There is significant over- 
lap between the distributions of velocity dispersions for 
the classical and pseudobulges. Nonetheless, it is clear 
that, in our sample, pseudobulges have on average lower 
velocity dispersions. We find in our sample that classi- 
cal bulges become exceedingly rare below central velocity 



dispersions of 100 kms ^. However, we caution that our 
sample is not volume limited. 

5.4. Velocity dispersion gradients 

Inspection of the individual rotation curves reveals a 
wide variety of structures; however, in particular the 
shape of the velocity dispersion profile seems to fall into 
two rough classes. In Fig. [6] we show the kinematic pro- 
files for the two galaxies NGIC 3898 and NGC 4448 from 
our sample. Depicted are the velocity, the velocity dis- 
persion, as well as the and moments of the Gauss- 
Hcrmite expansion of the LOSVDs. Dashed lines indi- 
cate the bulge radius from the photometric decomposi- 
tion. While in the case of NGC 3898 the velocity dis- 
persion rises all the way to the centre, NGC 4448 has a 
relatively flat dispers i on pr ofile within the bulge radius. 

Similarly to Fisher ( 19971 we examine, the logarithmic 
slope of the velocity dispersion within the bulge radius 
and call it 7. We derive the slope point-wise and then 
take the average, i.e. 



7 =< 



Alog(a) 
Alog(r) 



> 



T,<r<rb I 



(5) 



where rmin always excludes the inner FWHM of the see- 
ing of the particular observation and in some cases is 
chosen larger to exclude central features like nuclear re- 
gions of enhanced star formation (see Appendix pi). In 
Fig. [6] we also overplot lines which correspond to tne de- 
rived7 values. Further, in order to avoid a dependence of 
the slope on the particular binning scheme of each kine- 
matic dataset, we use a different binning for the purpose 
of determining 7: we bin radially in 5 equally-sized bins 
in log{r). In cases where the resulting bins do not all at 
least contain one data point, we use our previous bins. 

An alternative to the presented method is using the ra- 
tio of the averaged velocity dispersions within two annuli 
within the bulge radius 



< ^ >|ri,/3<r<r6 



(6) 



as proxy for the slope. The choice of r^/S as cut radius 
for the two different annuli is somewhat arbitrary, but 
we do not find a strong dependence of our results on the 
specific radius chosen. Both values for the slope, 7 and 
5, are reported in Tab. [9] 

We find that all bulges which are classified as pseudo- 
bulges indeed show fiattened velocity dispersion profiles 
or even sigma droro (e.g. NGC 3351, NGC 3368, and 
NGC 3627 in Fig. 16). The dispersion profiles of many 



pseudobulges are sometimes slightly asymmetric. On the 
other, hand a majority of the classical bulges show cen- 
trally peaked velocity dispersion profiles (e.g. NGC 1023, 
NGC 2841, NGC 2880, and NGC 3245 in Fig.[l6]). 

Fig. [7] summarizes this finding qualitatively, where we 
plot the velocity dispersion profiles along the major and 
minor axes for all our bulges separated by bulge type, 
normalized by central dispersion and bulge radius. We 
do not plot bulges that were left unclassified. For this 
plot we adjust the bulge radius that was obtained from a 
major axis profile by the mean ellipticity in the bulge re- 
gion. Whilst not as clear, partly due to the lower number 
of profiles, but partly probably also due to the subtleties 
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TABLE 8 

Format example of the measured stellar kinematics. 



Galaxy 


PA 


r 


V 


a 


/13 


hi 




[dog] 


n 


[kms^^] 


[kms^^] 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


NGC 1023 


87 


51.46 


186.55 ± 2.71 


92.16 ± 2.55 


-0.067 ± 0.017 


-0.037 ± 0.015 



Notes: 2J Identifier. 2) Observed position angle. 3) Distance from the center (positive: east; 
negative: west). 4) Velocity relative to systemic velocity. 5) Velocity dispersion. 6) Gauss- 
Hermite moment. 7) ^4 moment. 

The full listing is available electronicallyhttp : //cds .u-strasbg . f r/. 
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Fig. 5. — Histograms of tlie central velocity dispersions. The left panel discriminates bulge types by morphology, the right panel 
discriminates them by their Sersic index. 
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Fig. 6. — Major axis kinematic profiles for NGC 3898 and NGC 4448. Positive radii are east of the galaxy center. We plot from bottom 
to top the rotational velocity, velocity dispersion, and /14 moments. The curvature of the red lines correspond to the derived logarithmic 
slope of the dispersion profile, they are scaled to match the depicted profile. Their extend indicates the radial range which is taken into 
account for the derivation of the slope (see text) . 
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Fig. 7. — Major (top panels) and minor axis (bottom panels) ve- 
locity dispersion profiles, normalized by central velocity dispersion 
and bulge radius. Profiles of classical bulges are plotted in red in 
the respective upper panels, those of pseudobulges in blue in the 
corresponding lower panels. Different color shades correspond to 
different galaxies. The thick black line shows the median profile for 
all bulges in one panel. The bulge radii for the minor axis profiles 
have been corrected using the mean bulge ellipticity according to 
fb = (!-<£ >)-n. 

of choosing a correct radius for the normalisation, we 
again find that classical bulges tend to show centrally 
rising velocity dispersions. 

In Fig. |8] we now plot the Sersic index from the pho- 
tometric decomposition as a function of both metrics for 
the slope of the velocity dispersion. Similar to the distri- 
butions of central velocity dispersion, there is significant 
overlap in profile slope. Nonetheless, the bulges with 
large values of Sersic index tend to have steeply decay- 
ing dispersions profiles. Similarly the bulges with low 
Sersic indices more commonly have flat dispersion pro- 
files. This result is true for both the logarithmic slope of 
dispersion, and the dispersion ratio. 

The increasing slope of velocity dispersi on with Sersic 
index is not fully unexpected. For instance, Ciotti ( 1991 1 



describes a series of models for isotropic and spherical 
galaxies which have a surface brightness profiles that fol- 
low a Sersic law. He gives projected velocity dispersion 
profiles for his models and outside of the very central re- 
gions (r > Te/lO), and for Sersic indices larger than one, 
the slope is a monotonically increasing function of n. We 
calculate slopes and a ratios for these profiles in a similar 
manner as we did for our data. One caveat of this exercise 
is that our definition of a bulge radius is not applicable 
in the case of the one-component models. Also, we have 
to chose an inner cut radius for the fit as the models fea- 
ture central sigma drops in the case of small n. Sigma 



drops are a n observed phenomenon (e.g. Falcon-Barroso 
et al. 2006 1, but our spatial resolution is typically not hne 
enough to resolve these. We somewhat arbitrarily fit for 
7 in the radial range of Te/lO < r < and calculate cr 
ratios for rg/lO < r < r^/S and re/3 < r < r^. Note that 
our effective bulge radii are on average 15% smaller than 
the definition of the bulge radius that we use through- 
out the work. In Fig. [8] we overplot the obtained values 
as a black line. The dashed lines show the range of val- 
ues one would obtain by choosing 50% larger or smaller 
outer cut radii for the integration. While the spherical 
and isotropic galaxies are a very simplistic model for the 
variety of bulges in our sample, one can see that the 
general trends are reproduced, however a more detailed 
dynamical modelling is needed to confirm this result. 

In the major axis plots all unbarred pseudobulges fall 
below or very close to 7 = —0.05 and 6 — 1.06 (the cor- 
responding values of the isotropic models for n — 2) and 
only one unbarred classical bulge falls below 7 = —0.05. 
However 3 out of 8 unbarred classical bulges do fall signif- 
icantly below 5 = 1.09 suggesting that 7 is more success- 
ful in discriminating bulge types. Again this picture is 
complicated further once barred galaxies are taken into 
account. The additional component of a bar seems to 
lead towards flatter dispersion profiles. 

5.5. Influence of seeing on velocity dispersion 

The seeing disk and the width of the slit will smear 
the observed velocities and can create increases in the 
observed line of sight velocity dispersion. This effect is 
commonly known as slit smearing. All data presented 
here were observed with a slit width of 1 arcsecond. 
The effect of slit smearing on the velocity dispersion is 
therefore expected to be negligible compared to the ef- 
fect caused by the seeing (> 1.2 arcseconds in all cases). 
At least two galaxies NGC3384 and NGC3521 do show 
pea ks in velocity dispersion in the central arcseconds (see 
Fig. 16 1 . In both galaxies the velocity profile also rises 
rapidly in the centre. We test whether this rapid rise in 
combination with the seeing may be responsible for the 
observed dispersion peak. We model the point spread 
function (PSF) with a Gaussian of the same FWHM. 
We then calculate the standard deviation of the veloc- 
ity which is weighted by the PSF amplitude at all radii 
and subtracted the result from the observed velocity dis- 
persion. In this simple one-dimensional model the PSF 
smearing does generate a central peak which is of similar 
size and amplitude as the observed one. We cannot rule 
out the possibility that the central peaks of NGC 3384 
and NGC 3521 can be explained through PSF-smearing 
alone. We however refrain from correcting the presented 
velocity dispersions as an accurate correction has to in- 
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Fig. 8. — Sersic index n is shown as a function of both metrics for the flatness of the velocity dispersion profile, 7 and 5, respectively. 
Upper left: Shows the major axis logarithmic slope of the velocity dispersion. Red squares and blue circles correspond to classical and 
pseudobulges respectively. Open symbols label barred g alaxies. The green diamonds represent unclassified objects. The black solid line 
shows the respective behaviour of the isotropic models in |Ciotti| | |1991| , here the slopes were calculated in the radial range Ve/lO < r < re- 
The dashed lines show the 7 values one would obtain by choosing the outer integration radius 50% smaller or larger. The horizontal line 
marks a Sersic index of two. The vertical dotted lines mark 7 = —0.046 and 5 = 1.06, the respective values that the isotropic models take 
for a Sersic index of 2. Upper right: Shows the ratio of the averaged velocity dispersion in two different annuli. NGC3593 falls far to 
the left with 5 = 0.7 (see discussion in jj6]and Appendix [b|| . Lower panels: Same for the minor axis dispersions. The radii were adjusted 
according the mean bulge ellipticity fj, = (1— < e >) ■ r^. the error bars correspond to the typical errors in the derived quantities, they 
also apply to the upper panels. 

elude the knowledge of a high resolution luminosity pro- 
file and a more rigorous, 2-dimensional modelling of the 
PSF. We rather exclude the central peaks from the fur- 
ther analysis. 



5.6. Distribution of /13 and /14 moments 

As /i3 measures the asymmetric deviation from a 
purely Gaussian distribution it detects lower velocity 
tails of the velocity distribution al ong the lin e of sight. 
Such tails arise natu rally in disks ( Binney fc Tremaine, 
1987[ ) . [Bender et al.| ( |1994[ ) found that local /13 and local 
v/a are strongly anti-correlated wi th a slope of — 0.12 in 
their sample of elliptical galaxies. Fisher (1997) finds a 
similar anti-correlation in the inner regions of his lenticu- 
lar galaxies but also sees that, for a number of his objects, 



at values of w/cr « 1 the anti-correlation turns, at least 
briefly but abruptly, into a correlation. 

We reprod uce the plot for the local correlation of hs 
and v/a from'Bender et al. ( 1994 1 for our sample in Fig. M 
and color-code pseudobulges in blue and classical bulges 
in red. We find that the same correlation is reproduced in 
our intermediate type galaxies. The moments are gen- 
erally anti-correlated with v/a out to v/a Pi 0.5, irrespec- 
tive of bulge type. A linear fit to the complete set of data 
points gives a slope of (/13 = (—0.106±0. 001) -ti/cr). Sep- 
arate fits to the subsample of classical bulges and pseudo- 
bulges give values that are indistin guishable within the 



errors. A Kolmogorov-Smirnov test (Smirnov 1939 Press 

within Individ 



2002 1 for the median values oiv-a ^ -h^ 
ual galaxies yields a probability of 76% for the hypothe- 



Kinematics across bulge types 



15 



0.3 
0.2 
0.1 
-c" 0.0 
-0.1 
-0.2 
-0.3 




classical 
pseudo 
fit to classicals ' •■ 
fit to pseudos 
Bender et. al 94 for Es 



-2 



-1 





vja 



in bu l ges w ith larger averaged vja (Fig. 10). Bender 



et al. ( 1994 ) also discuss the possibility of a similar trend 
in their subsample of rotationally flattened galaxies. 

There is a mild indication that pseudobulges and clas- 
sical bulges show different distributions in the average 

< /i3 > and < /14 > moments. For pseudobulges we find 
an average value for < /13 > of 0.06 with an RMS scat- 
ter of 0.03 while for classical bulges the mean value of 

< /i3 > is 0.04 with a scatter or 0.03. A KS test and a 
Student's two-tailed t-test for two independent samples 
yield a probability of 0.3% and 3%, respectively, for the 
two subsamples to stem from the same distribution. For 

< /14 > we find and average of 0.04 with a scatter of 0.02 
in the pseudobulges, and 0.02 with a scatter of 0.02 in 
the classical bulges. Here the KS test finds a 5% prob- 
ability for the null hypothesis while the t-test yields a 
1.5% probability. As both the ^,3 and the moment are 
affected by the inclination and the scatter is large, this 
trend has to be taken with caution. A larger sample and 
kinematic modelling will be needed to confirm if this is a 
signature of systematic different anisotropics in the two 
classes of bulges. 



Fig. 9. — Local correlation between h-j, and v/cr along the major 
axis for the galaxies in our sample. We only plot points for which 
the error in /13 is lower than 0.05. Plotted are all galaxies for which 
the bulge was classified either as classical (red) or a pseudobulge 
(blue). Typical error bars are shown in the upper right of the 
diagram. The red and blue lines correspond to the fitted linear 
correlations for the classical and pseudobulges respectively. The 
dashed black liri e represents the value for the correlation that |Ben^ 
|der et ah] | |1994[ l obtained for their sample of early types. 



sis that the classical and the pseudobulges stem from the 
same distribution. This local correlation is rep roduced 
in the mean values for the bulge region (Fig. 10 1. 

We further test for a possible correlation with _ff-band 
bulge magnitude (see Fig. 10) and the bulge averaged 
value < /13 >. 



We do not see any correlation between 
bulge luminosity and < hs >. 

The hi moment of the Gauss-Hermite expansion mea- 
sures the symmetric deviation from a Gaussian distri- 
bution. Negative /14 describe a more boxy, centrally 
flattened distribution, more positive values describe cen- 
trally peaked distributions with extended wings. The 
averaged ft.4 moments in the bulges are generally close 
to zero, the median for the complete sample of major 
axis spectra is 0.03 with a standard deviation of 0.046. 
None of our bulges sh ow obvious dips in the /14 profile as 
the ones described by Debattista et al. (2005); |Mendez- 



Abreu et al. ( 2008 ) . But this diagnostic for boxy-peanut 
shape bulges only applies to low inclinations {i < 30°), 
given that the inclinations of most of our galaxies is 
larger than 30° (41 out of 45) this is not further sur- 
prising. However, fourteen galaxies show a double peak 
in the /14 profile within the bulge region (e.g. NGC 1023, 
NGC3031, NGC 3945, and NGC 7331 in Fig.ITel) . This 



is typically seen in combination with a rapid increase of 
the rotational velocity and relatively strong m omen ts. 
From our Monte Carlo simulations described in §4.1| we 
can rule out that the observed peaks are a result of a 
degeneracy between /13 and /14 moments in the fit. 

We find no correlation between the averaged < /14 > 
moments and the bulge luminosities. However, while the 
error bars are large, larger moments seem to be found 



5.7. Extreme moments and multiple kinematic 
components 

Five galaxies show extreme /13 and /14 moments. The 
most extreme case, NGC 3521, (a classical bulge) exhibits 
values of /13 and /14 as large as 0.24 and 0.35, respec- 
tively (see Fig. [16]). NGC 3945, NGC 4736, NGC 7217 
(all pseudobulges)7 and NGC 7331 (unclassified ) show 
values of /13 and /14 of up to 0.2. The LOSVD is poorly 
reproduced by a Gauss-Hermite expansion at t he r adii of 
such extreme higher moment values (see Fig. 11). The 
reason lies in the existence of a secondary kinematic com- 
ponent in al l those cases. F o r NG C 3521 this has been 
reported by Zeihnger et al. (2001) who attributed the 
counter-rotating stellar component to the presence of a 
bar. The two-com ponent nature of NGC 7217 was dis- 
covered before by [Merrifield fc Kuijken ( 1994 ) . They 
suggest that the second component is the result of an 
extended period of accretion with intermitte nt change of 
angula r momentum of the infalling material. Prada et al. 
( 1996 ) reported a counter-rotating bulge in NG C 7331. 
I'wo mo re systems in our sample, NGC 2841 (Bertola 



fc Corsini|il999) and NGC 3593 ( [Bertola et arfTUM^ , 
were reported to host counter-rotating components, and 



NGC 4698 



Falcon-Barroso et al 



( Corsini et al 



1999 


Bertola et al. 


1999 


see P 


izzelia et al. '(2004) tor 



which the rotation curve itself already reveals counter 
rotation through a twist — the rotation changes sign 
with respect to the systemic velocity at a radius of about 
20 arcseconds (see Fig. 16 and Appendix [b|) . The case is 
similar whilst not as pronounced for NGC 4698 where 
the rotation _curve becomes very flat towards the cen- 
ter (see Fig. 



16) 



However, in these cases the secondary 
component does not result in unusually strong /13 and 
moments in our data. 

In an attempt of a fairer treatment of their complex- 
ity we decomposed the FCQ-derived LOSVDs into two 
separate Gaussian components in a similar manner to 
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Fig. 10. — Major axis correlations between bulge averaged Gauss-Hermite moments < /13 > and < /14 >, bulge luminosities and < v/a >. 
Upper left: < h-j > a,s function of bulge ff-band magnitude. Red circles represent classical bulges, blue circles are pseudobulges, green 
circles represent unclassified bulges. Upper right: < /i3 > as function of < v/cr >. No inclination corrections were applied. Lower left: 
< ?t4 > as function of bulge magnitude. Lower right: < ?t4 > as function oi < v/a >. Representative error bars are displayed in the lower 
right hand part of each panel. In the case of the bulge magnitude they correspond to the typical error. For all other plotted quantities the 
larger error bars correspond to the typical RMS scatter of that quantity within the considered radial range, whereas the smaller error bar 
corresponds to the formal error of the derived average. 



We 



Scorza fc Bender| ( |1995[ ); |Zeilmger e tlT (2001') 

the Metrop olis- riastings algoritiim ( Hastings_1970j Press 



et al.][2007 ) to infer the parameters and error bars. Be- 
tore tne ht, the spectra are binned to a minimum S/N of 
75 per pixel. We run four simultaneous chains for each 
radial bin. The step width is tuned to achieve a 25% 
acceptance ratio and after convergence the first half of 
the chain is discarded (clipped). The run is aborted if 
the chains do not converge after 100.000 steps. 

The 
and 



Fig. 



12 shows the result of this decomposition, 
plotted values are the maximum-likelihood values, 
the end of the error bars mark the 20% and 80% quan- 
tiles in all four chains after clipping. Central values with 
strong degeneracies between the parameter sets are omit- 
ted. In all five galaxies we do find significant second 
components under the assumption th at in dividual com- 
ponents are purely Gaussian. In Tab. [TO] we list the in- 
tegrated fractions of hght in the two dmerent kinematic 
components and compare those to the values that one 
would expect from the photometric decomposition. 



6. DISCUSSION 

6.1. Dichotomous dispersion profiles of classical and 
pseudobulges 

It is commonly assumed that the bulge light — typi- 
cally determined from a bulge-disk decomposition — rep- 
resents a dynamically hot component. Yet, it has been 
kno wn for a long time that rotation-supported bulges ex- 
ist ( [Kor mendy 1982). Also, many bulges have lower cen- 
tral vel ocity dispersions than expected fro m the Fabcr ^ 



Jackson, (,1976^) relatio n ( Kormendy fc Kennic utt„ 2004 [) 
I'alcon-Barroso et al.l ( |2006| ); IGanda et al.| ( |2006p hnd 
that several galaxies m the SAUKUJN survey have cen- 
ters that are dynamically colder than the surrounding 
disks. Thus, it is now clear that not all galaxies fit the 
picture that a bulge is a dynamically hot component. 

There is an observed dichotomy in bulge properties 
including Sersic index, bulge morpholog y, star forma- 



tion fc ISM properties, and optical color ( Carollo et al. 
1997"; 'Gadotti fc dos Anios' |200Tl [Ko rmendy Kenm" 
cutt 2004; Fisher 2006; Fisher fc D rOTy 2008 2010^"~Ri 



cently, F^iier & Drory ^ 2010[ ) showed that tlie dichoto- 
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TABLE 9 

Structural and kinematic parameters. 



Galaxy '■^b=fjrf < tb > < Cd > o"re/lo 7A/.7 Smj jmn Smn ^^iy 

arcsec arcscc kms^^ 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 



NGC 1023 


21.3 


19 








22 





55 


212.9 ± 5.2 


-0 


09 


1 


12 


-0.17 


1.176 


0.165 


NGC 2460 


8.3 


6 


6 





19 





25 


111.4 ± 3.5 


-0 


01 





99 


0.10 


1.026 


0.358 


NGC 2681 


14.6 


13 


2 





11 





20 


112.5 ± 1.3 


-0 


16 


1 


19 








NGC 2775 


19.4 


16 


8 





10 





16 


173.9 ± 13.7 


-0 


04 


1 


02 


1.03 


1.027 


0.133 


NGC 2841 


17.4 


15 


2 





22 





49 


222.2 ± 19.3 


-0 


11 


1 


15 


-0.09 


1.124 


0.207 


NGC 2859 


30.0 


27 


6 





16 





22 


176.8 ± 5.4 


-0 


06 


1 


11 


-2.04 


1.379 


0.678 


NGC 2880 


26.7 


22 


7 





20 





35 


142.2 lb 5.3 


-0 


21 


1 


35 


-0.22 


1.394 


0.496 


NGC2964'' 


3.4 


3 


1 





15 





28 


88.4 ± 1.3 














0.117 


NGC 3031 


72.0 


61 


3 





24 





43 


157.5 ± 13.6 


-0 


06 


1 


03 


-0.03 


1.006 


0.400 


NGC 3166 


9.9 


9 


1 





39 





25 


151.4 ± 6.1 





00 


1 


04 


0.05 


0.981 


1.365 


NGC 3245 


9.5 


8 


5 





20 





44 


225.2 ± 8.3 


-0 


21 


1 


27 


-0.22 


1.184 


0.300 


NGC 3351 


14.2 


12 


9 





16 





24 


90.0 ± 4.2 





05 


1 


01 


-0.05 


1.015 


0.778 


NGC 3368 


23.3 


20 


4 





17 





34 


122.5 ± 6.6 


-0 


01 





99 


0.03 


0.883 


0.574 


NGC 3384 


15.5 


14 


4 





20 





34 


150.3 ± 2.4 


-0 


09 


1 


10 


-0.13 


1.150 


0.521 


NGC 3521 


12.2 


10 


8 





35 





45 


129.5 ± 2.9 


-0 


01 





99 


-0.14 


1.098 


0.905 


NGC 3593 


32.1 


29 


5 





49 





62 


62.3 ± 3.1 





20 





70 






0.065 


NGC 3627 


11.7 


10 


9 





27 





51 


116.1 ± 3.9 


-0 


05 


1 


03 


0.01 


1.003 


0.267 


NGC 3675 


9.8 


8 


5 





29 





49 


114.7 ± 5.3 


-0 


03 


1 


03 






0.363 


NGC 3898 


17.7 


15 
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25 





41 


219.0 ± 8.3 


-0 


14 


1 


21 






0.295 


NGC 3945 


33.3 


31 








19 





17 


183.1 ± 5.4 


-0 


03 


1 


12 


-0.11 


1.088 


1.062 


NGC 3953 


16.6 


14 


6 





26 





48 


110.6 ± 3.1 





02 





95 






0.178 


NGC 3992 


14.9 


13 


2 





22 





49 


144.2 ± 9.5 





07 





97 






0.406 


NGC 4030'' 


3.6 
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11 





19 


102.9 ± 4.5 














0.213 


NGC 4203 


16.4 


14 


7 





11 





11 


170.1 ± 3.6 


-0 


18 


1 


31 






0.249 


NGC 4260 


9.0 


7 


3 





21 





53 


143.8 ± 14.3 





03 





98 






0.131 


NGC 4274 
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11 


3 





40 





34 


106.9 ± 5.3 





08 





90 






0.842 


NGC 4314 
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8 


6 
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45 


123.3 ± 5.1 


-0 


01 


1 


00 






0.825 


NGC 4371 


25.5 


22 


9 





29 





33 


125.8 ± 5.0 


-0 


12 


1 


20 


0.00 


0.983 


0.426 


NGC 4379 


10.0 


8 


6 





11 





20 


121.0 ± 4.6 


-0 


20 


1 


16 






0.183 


NGC 4394 


15.2 


14 


1 





12 





37 


80.0 ± 3.1 


-0 


01 





91 






1.451 


NGC 4448 


9.6 


8 


5 





26 





43 


98.5 ± 3.7 


-0 


03 


1 


03 






0.446 


NGC 4501 


7.0 


6 


2 





19 





45 


144.2 ± 4.9 





00 


1 


01 


0.07 


0.937 


0.390 


NGC 4536 


10.9 


10 


1 





39 
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° Structural and kinematic parameters for the galaxies in our sample. 1) Target name. 2) Radius of 
equal bulge and disk surface brightness. 3) Adopted bulge radius for this study. 4) Mean apparent bulge 
elliptieity. 5) Mean apparent disk ellipticity. 6) Central velocity dispersion averaged within one-tenth of the 
bulge effective radius. 7) Slope of major axis velocity dispersion profile. 8) Major axis ratio of the averaged 
velocity dispersion within the annulus r < r^/S to averaged dispersion within rf,/3 < r < r(j. 9) Slope of 
minor axis ve locity dispersion profile. 10) Minor axis ratio of velocity dispersion. 11) < •u'^ > / < cr^ > 

! Binney|2Q05[ for the inclination corrected velocity. 
comments: a) The low inclination of this galaxy prevents us from deriving an inclination corrected velocity 
and hence < > / < a'^ > . b) The bulge is not sufficiently resolved to calculate the slopes of the velocity 
dispersion, c) Surface brightness does follow a typical bulge/disk profile. We do not decompose the profile 
and only present the kinematics data here, d) No decomposition, this is the innermost value. 



mous properties in Sersic index, morphology, and ISM 
properties are consistent. Furthermore they show that 
bulges of different type occupy different regions in the 
projection of fundamental plane properties, thus indi- 
cating that there are very likely two physically distinct 
classes of bulges. Are dynamics part of this dichotomy? 

The high spectral resolution of 39 kms~^ of our data, 
enables us to recover dispersions out into the disk regions 
in many of our targets — a feature uncommon to many 
similar surveys. We extract LOSVDs using the FCQ al- 
gorithm with additional procedures to account for neb- 
ular emission and template mismatch. We recover v, a, 



/13, and /14 moments of a Gauss- Hermite model of the 
LOSVDs as function of radius. 
We observe a gr eat variety of shapes of kinematic pro- 



files (see Fig, 

pendix^. S i rnTIar ly to |Falc6n-Barroso et al 



16 and their detailed descript i on in Ap- 
- - ■ - - p006| and 



Uanda et'al. ( 2006 ) , we hnd that it is not necessarily true 
that the center of a bulge has the highest observed veloc- 
ity dispersion (e.g. NGC 3593). In our sample only ^1/3 
of the galaxies have centrally peaked velocity dispersion 
profiles (like NGC 3898). Many galaxies have roughly flat 
velocity dispersion profiles. In these galaxies there is no 
apparent transition in velocity dispersion from the bulge 
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Fig. 11. — Example for the double-Gaussian decompositions for 
NGC3521. The FCQ derived full line of sight velocity distribu- 
tion for five selected radii is plotted in black. The two-Gaussian 
kinematic decompositions plotted in red and blue. 



TABLE 10 

Light fractions in kinematic subcomponent. 



galaxy 


morph. 


disk light 


light in fast component 






[%] 


[%] 


(1) 


(2) 


(3) 


(4) 


NGC 3521 


SAB(rs)bc 


34 


44 


NGC 3945 


(R)SB+(rs) 


2 


34 


NGC 4736 


(R)SA(r)ab 


55 


50 


NGC 7217 


(R)SA(r)ab 


55 


38 


NGC 7331 


SA(s)b 


26 


47 



no es. Comparison of the amount of light in the lower dispersion 
component to the light that one would expect from the photomet- 
ric decomposition. We integrate the amplitudes of the G aus sian 
fits to the two components over all radii shown in Fig. |12| 1) 
Galaxy. 2) Morphology (RC3). 3) Integrated light in the disk 
component from the extrapolation of the outer disk exponential 
profile. 4) Integrated light in the lower dispersion component. 



to the disk region unlike in the stellar surface brightness 
profile (e.g. NGC 4448 and NGC 5055). 

It is interesting to note that to the radial extent of our 
data, in many of these cases the disk velocity dispersion 
is as high as the central velocity dispersion of the galaxy. 
For example in NGC 4448 the disk velocity dispersion 
remains above 100 kms^^. Therefore some disks of spiral 
galaxies are not necessarily cold stellar systems over the 
radii that we cover in this study. 



From minor axis data of 19 SO to Sbc bulges (Falcon- 
Barroso et al.|2003 ) showed that higher ellipticity bulges 



have shallower velocity dispersion profiles profiles. If 
pseudobulges appear photometrically flattened then one 
might expect that the steepness of the dispersion pro- 
file should correlate with bulge type. We show in Fig. [7] 
that the shape of the velocity dispersion profile corre- 
lates very well with bulge type. Galaxies with classi- 
cal bulges have centrally peaked profiles. Galaxies with 
pseudobulges have, on average, flat dispersion profiles. 



We have attempted to quantify this using the logarith- 
mic derivative of the velocity dispersion as function of 
radius (Eq. Island also the ratio of dispersions at differ- 
ent radii (EqlOl . We find that pseudobulges and classical 
bulges occupydifferent regions in the parameter space of 
logarithmic derivative of velocity dispersion and Sersic 
index (see Fig. |8]) in a way that is not inconsistent with 



models of dynamically isotropic systems ( Ciotti|[r991[) . 

It is important to note that the dynamics of a few 
galaxies are not well described by a simple monotonic 
trend of velocity dispersion with radius; we stress that 
for the purpose of this paper we are interested in the bulk 
properties of the distribution of stellar dynamics. The 
great variety of shapes in dispersion profiles that we ob- 
serve (for a detailed description see Appendix |b| is likely 
to be a consequence of the fact that there are multiple 
ways to heat galactic disks, for example through merg- 
ers [van Albada,, 19821 IQ uinn et al.|ll993| lEliche-Moral| 
et al.|2006l iHopkms et al. 2008), bars (Saha ct al. 2010r 



and other disJt insta bilities ( jSellwood fc Wilkinson„1993"j 
Combes et al.||199"0 ). Even under the strong assuniption 
that classical and pseudobulges are dynamically distinct, 
it does not seem plausible that any simple description 
of the kinematic profile cleanly separates classical from 
pseudobulges. 

A few classical bulges in barred galaxies such as 
NGC 3992 do not seem to fit this general picture. How- 
ever, all those galaxies are barred and it is conceivable 
that bars may distort the kinematic profile of a classi- 
cal bulge as they vertically heat the disk they reside in 



(Gadotti & de Souza 2005 Saha et al. 2010). Central 



velocity dispersions lie higher by a factor of two than at 
the bulge radius in the most extreme cases in our sam- 
ple, e.g., in NGC 1023, NGC 3898, and NGC 4203. If bars 
raise velo city dispersions by a factor of up to four as sug- 



gested by Saha et al. (20101, the signature of a central 
dispersion peak can of course be easily washed out. 

NGC 4826 has extreme amounts of dust in its central 
region — hence also its name black eye or evil eye galaxy. 
The bulge was consequ ently classified as pseudobulge by 
[Fisher fc Drory (2008). It stands out however, as it has 
a relatively high Sersic index of 3.9 ± 0.88 for a pseudo- 
bulge. The V-band value of n = 3.94 ± 0.68 (iFisher & 



Drory|2008 ) agrees well. The bulge radius of 25.4 arcsec 
onds seems small once the kinematic data are taken into 
account. The velocity dispersion starts rising at about 
50 arcseconds already, which corresponds to the radius of 
the final fiattening of the rotation curve. It is conceivable 
that the large amount of dust in its center, which is easily 
also visible in the infrared, may affect the decomposition. 
If one was to take the value of 50 arcseconds as the bulge 
radius then the 7 value would become -0.33 (0.01 before) 
and the sigma ratio would take a value of 1.29 (1.08 be- 
fore) . This would move NGC 4826 significantly further to 
the right in both plots into the region occupied by clas- 
sical bulges, in much better agreement with the Sersic 
index. Further the disk of NGC 4826 is relatively free of 
dust and actually resembles an SO. We hypothesize that 
the unusual morphology is a result of a recent merging 
event. A satellite may have fallen into an SO-like disk 
and brought in dust and triggered star formation. This 
hypothesis is supported by the exis tence of two counter 



rotating gaseous disks observed by Braun et al. (1992 



19941. We labelled NGC 4826 as non-classitied througE 



Kinematics across bulge types 



19 



l.OOF 

0.75 

0.50 

0.25 

0.00b 



NGC3521 major axis 




-10 10 

r [arcsec] 



1.00 
0.75 
0.50 
0.25 
0.00 

250 



a 



100 

300 



NGC3945 major axis 





L..'i.i-,.,: : 








t r — i : 




r [arcsec] 



1.00 
0.75 
0.50 
0.25 
0.00 
200 
150 
100 

50 
200 
100 


-100 
-200 



NGC4736 major axis 



-40 




-20 20 
r [arcsec] 




40 



^ 200 
> 100 

a 0- 

"-100 T 

* -200 - 



NGC7217 major axis 




-40 -30 



-20 -10 10 20 30 40 
r [arcsec] 



1.00 
0.75 
0.50 
0.25 
0.00 

- 300 
a 200 

" 100 

b 



300 
200 
> 100 

^ -100 
s -200 
-300 



NGC7331 major axis 



; I 












r; 







-30 -20 -10 10 20 30 
r [arcsec] 



Fig. 1 2. — Kinematic decompositions of the major axis data of NGC3521, NGC3945, NGC4736, NGC 7217 and NGC7331. In a procedure 
similar to |Scorza &: Bender] | |l995[ ; [Zeilinger et al.| | |2001[ | we fit two Gaussians to the FCQ derived full LOSVDs at all radii. Blue and red error bars 
show mean velocity and dispersion of the two Gaussians where the Metropolis-Hastings algorithm chains achieved convergence. The length of the 
error bars represent the 20% and 80% quantiles. Black curves indicate the FCQ fitted moments of a single Gauss-Hermite expansion of the LOSVD. 
Here the error bars are comparable to the thickness of the line. The upper panel shows the relative weight of the fast component with respect to 
the total light. The dashed curve shows the disk to total ratio from the photometric decomposition. Vertical lines indicate the bulge radius. 
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out this work. NGC 3593 has a very large value of 7 and 
a very small sigma ratio, i.e. it falls far to the left in 
both diagrams. This is a result of the strong depression 
in velocity dispersion in the bulge region. NGC 3593 is 
the only galaxy in the sample where we observe actual 
counter rotation in the sense of a change of sign of the 
mean rotational velocity in the bulge region. The small 
Sersic index of 0.81 supports the picture that the bulge 
region is dominated by a kinematically c old and distinct 
but luminous disk ( [Bertola et "aI1|1996[ ). NGC 2681 is 
classified as a pseudo bulge by morphology and yet has a 
relatively large Sersic index of n = 3.8. Further, it has a 
centrally peaked dispersion profile with 7 = —0.16 and 
S — 1.19. This agreement of photometric structure and 
dispersion slope prompts us to reassess the morpholog- 
ical classification. While the disk shows relatively little 
amounts of dust, a high contrast dust spiral within in 
the bulge easily seen in HST F555W, offers a clear sense 
of rotation. Also the spiral is not obviously misaligned 
with the outer disk. This galaxy may represent a pro- 
totypical case for the breakdown of the morphological 
classification sche me. However, it shows multiple bars 
— possibly three (Erwin & Sparke 1999) — and hence 
the central heating may also be a consequence of its com- 
plicated dynamical structure. NGC 3521 has a seemingly 
a relatively flat dispersion slope with values of 7 = —0.01 
and S = 0.99. We discuss this object at the end of the 
next Section. 

6.2. Rotational support 

In order to study the level of rotational support of 
a stellar system, it has become common practice to 
study its location in the VmaxI < u > vs. e dia- 



gram (Illingworth 1977 


Binney & Tremaine 


1987 


Ko- 


rmendy||1993|) — where 


"^max measures the maximum 



persion within a certain radius, and e the system's el- 
lipticity. One can also directly compare Vmax/ < c > 
to the expected values of an oblate-spheroidal system 
with isotropic velocity disp ersion. For instance, |Kor-| 
mendy & Illingworth ( |1982 ) define the anisotropy para- 



meter (v/a)* = {vmax/ < (J >)/ \J (-1 (1 — e) as a measure 
for the rotational support of a stellar system. Values of 
{v/a)* ss 1 point towards a support by rotation whereas 
values < 1 indicate support by anisotropy. Those mea- 
sures involve the ellipticity of the system which is typi- 
cally subject to relatively large uncertainties, especially 
when measured for galaxies which are dominated by large 
quantities of dust. Here, we decide to rather examine 
the local {vcorr / i-e. as a function of radi us, and 
the a veraged values of < Vcorr^ > / < cr^ > (Binney 
20051 across the bulge region, Vcorr = Vobs/ sin(i) is the 
inclination corrected velocity at a given radius. We use 
inclinations from Hyperleda (see Tab. [T]). We apply no 
further co rrec tion to the velocity dispersion. 

In Fig. 13 we plot {vcorr/cr){r) separately for classi- 
cal and pseudobulges. We further plot histograms of the 
bulge-averaged quantities in Fig. |14| Again we normal- 
ize the radii by the bulge radius and exclude the cen- 
tral seeing FWHM from the analysis. While there is 
significant overlap between the two sub samples, pseudo- 
bulges are biased towards larger {vcorr / This is 
especially seen in the histograms for the averaged val- 
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Fig. 13. — Local {vcorr / cF)(r) along the major axis, the radii 
are normalized by bulge radius. The velocities are corrected for 
inclination through l/sin{i). Classical bulges are plotted in the 
upper panel, pseudobulges in the lower panel. The solid black line 
marks the median of all profiles. 



ues. A Kolmogorov-Smirnov test ( Smirnov|1939 Press 
2002| yields a probability of 0.8% (1.9%) tor the classical 



and the pseudobulges in the full (non-barred) sample to 
stem from the same distribution. A Student's two-tailed 
t-test for two independent samples yields a probability 
of 0.7% (2.5%) for the classical and the pseudobulges 
in the full (non-barred) sample. This result supports a 
picture of increased rotational s upport of pseudobulges 
that was originally described by Kormendy (1993) and 
disc ussed in detail in 'Korme ndy (S~kennicutt i2004 ), see 
also Kormendy fc Fisher^ ( ,2005 ) . 



> 



While the average values of 7 and < > / < cr^ 
are different for classical and pseudobulges, neither of 
the two quantities separates the bulge types. In Fig. 15 



we combine both and plot < > / < > against 
the logarithmic slope of the velocity dispersion, 7, for 
the bulges in our sample. In the left panel we discrim- 
inate bulges morphologically and in the right based on 
the bulge Sersic index. The dashed line is drawn to con- 
tain all the pseudobulges. Classical bulges both tend to 
have low < > / < ct^ > and steeper negative slopes 
than pseudobulges. 

There are several galaxies with large Sersic index that 
fall dynamically into the region of pseudobulges. We 
have discussed the possibility that bars mask or destroy 
the kinematic signature of a classical bulge. However two 
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Fig. 14. — Histog rams for the averaged values oi < > / < 
> |Binney||2005 I. The velocities are corrected for inclination 
through l/sinit). tilassical bulges are plotted in the upper panels, 
pseudobulgcs in the lower panels. 

non-barred high-Sersic index galaxies remain (NGC 3521 
and NGC 7331; NGC 3521 also has classical bulge mor- 
phology). Both galaxies, have centrally peaking veloc- 
ity dispersion profiles but they also show a rise in ve- 
locity dispersion at larger radius (see Appendix [A]) , yet 
still within the radius of the bulge. This is also a fea- 
ture that is prominently seen in the dispersion profile of 
NGC 3031, but in this case the effect on the position in 
the < > / < (7^ > vs. 7 plane is not as dramatic. This 
behaviour is not what we observe in pseudobulgcs, which 
have fiat dispersion profiles. None of these three galaxies 
violate the general dichotomy observed in Fig. [Tj There- 
fore it is likely that their outlier location is due to a fail- 
ure in the machinery of measuring dynamical quantities. 
NGC 7331 and NGC 3521 both show signs of counter ro- 
tating components and NGC 3031 is well known to be 
interacting. It is possible in each of these galaxies that 
an outside mechanism is superimposing extra kinematic 
structure that is visible as an outside rise in the velocity 
dispersion profile of the bulge. We remeasure 7 changing 
the outer radius cut isolating only the central dispersion 
peak. In the case of NGC 3031 this moves the measured 
dynamical quantities into the region of parameter space 
that is only occupied by classical bulges; NGC 3521 and 
NGC 7331 both move significantly closer. 

Fig. [15] illustrates agreement between kinematic diag- 
nostics of the bulge dichotomy with structural and mor- 



phological indicators of bulge types. 

6.3. Multiple kinematic components 

We observe counter rotation seen as secondary com- 
ponents in the full shape of the LOSVD in five sys- 
tems (NGC 3521, NGC 3945, NGC 4736, NGC 7217, and 
NGC 7331; see ^5.7). It is striking how clearly the 
LOSVDs can be decomposed into a low-dispersion and a 
high-dispersion component in all these cases. It is tempt- 
ing to interpret the latter as the bulge and the former 
as the disk. However, if we plot the local disk-to-total 
ratios as obtained from the photometric decomposition 
over the values obtained fro m t he kinematic decompo- 
sition (upper panels of Fig. [l2| ), then we see that the 
disk contribution from the photometry falls short in all 
cases. Therefore, the observed low-dispersion component 
within the bulge region is not simply the extension of the 
outer disk as more light contributes to this component 
as one would expect from the extrapolation of the outer 
disk exponential profile alone. It is important to point 
out that a Gauss-Hermite distribution with moderate ^13 
moments can be modelled rather well by two individual 
Gaussians. /t3 m oments, though, occur naturally in disks 
(Binney & Tremainc 1987) and are not a signature of an 
actual second^component. Only in cases where the sec- 
ond component is clearly seen as a second peak in the 
LOSVD it is safe to assume that actually two distinct 
components contribute. 

The observation of a counter rotating comp onent is 
very interes t ing in the light of the findings of 'ElicheJ 
Moral et al. (2011 1. They use coUisionless A^-body simu- 
lations to study the characteristics of inner components 
such as inner disks and inner rings that were formed 
though a minor merger. For this, they simulate a number 
of mergers with different mass ratios and orbits. In gen- 
eral, while all their mergers formed an inner component 
supported by rotation, none of their mergers produced 
a significant bar. In their simulations all mergers with 
satellites on retrograde orbits do form a counter rotating 
component while none of them leads to actual counter 
rotation in the sense of a change of sign of the mean ro- 
tational velocity. A central increase of velocity dispersion 
and strong /13 and moments are observed in a major- 
ity of their models. They further find, however, that 
V /a and /i3 are generally anti-correlated within the inner 
components throughout all of their simulations. We see 
all these features refiected in the aforementioned galax- 
ies. Of particular interest to us is the double peak in 
that their model M6PlRh produced. Such double peaks 
are very pronounced in NGC 3945 and NGC 4736, but 
are also visible in the case of NGC 3521 and NGC 7331. 
In our small sample 3 out of the 5 systems for which 
counter rotation is observed are not barred, we think 
that a merging event is a likely formation scenario. 



7. SUMMARY 

In this paper we present kinematic profiles for the ma- 
jor axis of 45 intermediate type (SO-Scd) galaxies. Our 
survey differs from other similar surveys in that we are 
able to resolve lower velocity dispersions which allows us 
to study kinematic features in cold systems like disks and 
pseudobulgcs. We combine these data with bulge-to-disk 
decompositions of the stellar light. 
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Fig. 15. — The bulge averaged quantity < > / < cr^ > as function of the slope of the velocity dispersion, 7. In the left panel we 
discriminate bulge types based on morphology. Pseudobulges are colored blue, classical bulges red, and bulges that we do not classify are 
colored green. Light-shaded symbols represent barred galaxies or galaxies hosting an oval. In the right panel we discriminate by Sersic 
index (blue n < 2.1, red otherwise). For NGC3521, NGC3031, and NGC 7331 we remeasure 7 changing the outer radius cut to isolate 
only the central dispersion peak as described in the text. The corresponding new locations are marked as black open symbols. The dashed 
line is drawn to contain all the pseudobulges. 

We find that bulges that have increased rotational sup- 
port, as measured by larger values oi < v"^ > / < >, 
are likely to have lower Sersic indices and show disk-like 
morphology. 

Classical bulges on average tend to have higher central 
velocity dispersions than pseudobulges. In our sample 
the lowest central velocity dispersion in a galaxy with 
evidence for a classical bulge through a Sersic index of 
3.7 is arjiQ = 85 ± 2 kms'^ (NGC 7743). 

We observe — for the first time — a systematic agree- 
ment between the shape of the velocity dispersion pro- 
file and the bulge type as indicated by the Sersic index. 
Classical bulges have centrally peaked velocity dispersion 
profiles while pseudobulges in general have flat dispersion 
profiles and even at times show drops in the central ve- 
locity dispersion. We confirm that this correlation holds 
true if visual morphology is used for the bulge classifica- 
tion instead of the Sersic index, as it is expected from the 
good correlation between bu lge morphology and Sersic 



bulge-averaged quantities. We observe no correlation of 
the higher moments with bulge luminosity, however we 
find a weak correlation between the average values of 
and v/a. 

Through examination of the figures in Appendix [A] it 
is clear that the kinematic profiles of bulge-disk galax- 
ies commonly contain substructure. Furthermore non- 
axisymmetric features in the stellar structure such a 
bars make understanding the kinematics of these galaxies 
more difficult. Future progress will require 2D-methods 
capable of resolving low velocity dispersions commonly 
found in pseudobulges. We are currently e xecuting such a 
surve y using the VIRUS-W spectrograph ( Fabricius et al. 
2008|. 



index (Fisher & Drory 



2010) 



We observe that the disk regions of some of our galaxies 
have not always a low velocity dispersion. In some galax- 
ies the velocity dispersion remains above 100 kms"-'^ well 
into the region where the disk dominates the light. 

We confirm the previously described mu lticompo- 
nent nature of the full LOSVD in NGC 3521 peilinger 



[erar][200T| ) NGC 7217 ([M erri field fc Kuijken||l\)ij4| and 
NGC 7331 ^Prada et aT.||1996p and find two additional 
systems — namely IN GC 3945 and NGC 4736 - with sig- 
natures of multiple kinematic components. They be- 
come apparent through a secondary peak or pronounced 
shoulder in the full line of sight velocity distributions 
(LOSVD). We present double-Gaussian decompositions 
which show the presence of a counter rotating stellar 

component in all these sys tems. 

As in elliptical galaxies (Bender et al. 111994 1, we find a 



correlation of /13 and v/a, both locally as well as in the 
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APPENDIX 
KINEMATIC PROFILES 
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continued -Major and minor axis kinematic profiles for NGC 2859. We plot the SAURON results (de Lorenzo-Caceres 
in green. 
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continued -Major and minor axis kinematic profiles for NGC 2880. 
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(2001 1 were taken with a slit position angle of 157° whereas we observed at a position angle of 137° 
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continued -Major and minor axis kinematic profiles for NGC 3166. 
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continued -Major and minor axis kinematic profiles for NGC3368. We plot data from IHeraudeau et al.l (|1999|) in 
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continued -Major and minor axis kinematic profiles for NGC3521, see also figures [11] anci [12| 
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Fig. 16. — continued -Major axis kinematic profile for NGC3593, we plot the results of Bertola et al. (19961 in green. 
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continued -Major axis kinematic profile for NGC 3953. 
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Fig. 16. — continued -Major axis kinematic profile for NGC 3992. 
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Fig. 16. — continued -Major axis kinematic profiles for NGC4030. We plot the SAURON results of IGanda et all p006l 
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Fig. 16. — continued -Major axis kinematic profile for NGC4260. 
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Fig. 16. — continued -Major axis kinematic profile for NGC4274. We plot the SAURON results ( Falcon-Barroso et al.|2006 l in 
green. 
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Fig. 16. — continued -Major axis kinematic profile for NGC4314. We plot the SAURON results ( Falcon-Barroso et al.|2006 l in 
green. 
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Fig. 16. — continued -Major and minor axis kinematic profiles for NGC 7217, see also Fig. [T2| 
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continued -Major axis kinematic profile for NGC7331, see also Fig. |12| We plot results from [Vega Beltran et al^ 



(20011 for NGC7331 in green. 
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NOTES ON INDIVIDUAL GALAXIES 
Classical bulges 

NGC 1023 — .LBT-.. — Clean classical morphology in HST F555W. The rotation curve steps rapidly from w 
+60 kms""'^ to « -60 kms~^ in the central two arcseconds and rises then gradually to a value of sa 200 kms^^ at 
50 arcseconds. The pr ofile shows an eq ually fast change within the inner « 2 arcseconds from -0.03 to 0.03 with 
opposite sign. Like |Emsellem et al. ( 2004 ) we see that the v //13 anti-correlation turns into a correlation outwards of 



ilO arcseconds. The velocity dispersion profile rises all the way to the centre. Outside of 50 arcseconds the velocity 
dispersion profile flattens out at about fOO kms""'^ which coincides with a flattening in the rotation curve. This is 
significantly beyond the bulge radius of « 19 arcseconds. The minor axis rotation is mostly close to zero at larger 
radii but becomes negative inwards of 4 arcseconds -25 kms~^ at the centre). The acquisition image does not 
show an obvious offset of the minor axis slit but we note that due to the rapid rise of major axis rotation in the 
central arcseconds already a small offset of {~ 0.5 arcsecond) to the west suffices to explain the observed behaviour. 
The moments become positive in the same radial range, which is expected if the velocity offset is due to actual 
rotation. The minor axis moments show a double peak at ± 7 arcseconds with maximum values of /14 ~ 0.04. The 
continuously centrally rising velocity dispersion of the major axis is reproduced on the minor axis. 

NGC2775 — .SAR2.. — The HST F606W image shows a very clear classical morphology, the F450W image shows 
very little amounts of dust in the central region. We see a depression in the velocity dispersion profile inwards of 
5 arcseconds in the major axis profile as well as in the minor axis profile which coincides with a steeper part in the 
rotation curve. Eskridge et al. (2002) describe a large, slightly elliptical bulge which contains a bright nuclear point 
source. The Spitzer MlL'b 24 /im image shows a resolved nuclear source of emission which clearly stands out from a 
region of low emission which again in size roughly corresponds to the bulge radius. This may hint at a cold nuclear 
component with active star formation which dominates the kinematics. We exclude the data inwards of 5 arcseconds 
from our analysis but note that this choice does significantly affect the position of NGC 2775 in the plane of Sersic 
index versus velocity dispersion slope. 

NGC 2841 — .SAR3*. — The HST F438 image shows a weak nuclear dust spiral that is misaligned with the outer 
disc. The larger scale bulge morphology is smooth and shows little sign of dust and no spiral pattern. The /13 profile 
is anti-correlated with the rotation curve inwards of 4 arcseconds but then changes sign and becomes correlated with 
the velocity until about 20 arcseconds. The velocity dispersion profile is centrally rising and may show a little shelf 
inwards of « 5 arcseconds. The minor axis rotation curve shows an offset within the bulge radius. 

NGC 2859 — RLBR-I-.. — Prominent outer ring galaxy. This galaxy has no close-to V^-band HST image available. 
The bulge morphology in HST F814W and in the acquisition images is generally smooth and classical with few weak 
dust lanes. Our HET long slit data do not cover the full bulge region (r^ = 27.6 arcsec). Still the rotation curve 
starts to flatten out at our outermost data point at r « 8 arcseconds. Within this region the velocity dispersion 
profile rises centrally from about 125 kms~^ at ±6 arcseconds to 175 kms~^ in the centre. The /13 moments are clearly 
anti-correlated with velocity. The /14 moments show indication of the double peak signature at r w ±5 arcseconds. 
The minor axis kinematic data only reach out to 6 arcseconds with little rotation along the minor axis (less than 
25 kms"^). The coverage of the dispersion profile in insufficient to judge whether the central rise seen along the major 
axis is reproduced. The minor axis /13 and moments are somewhat noisy but do not exceed values of 0.05 and show 
no significant trends. 



NGC 2880— .LB.-.. — The bulge morphology is classical in HST F555W. |Erwin| ([2004') finds indication of a weak 



inner disk but acknowledges that this is the weakest case in his sample. The velocity dispersion profile rises centrally 
with a weak non-symmetric shelf-like structure inwards of w 4 arcseconds. The /13 moments are anti-correlated with 
velocity within the bulge. 

NGC 3031 — .SAS2.. — M81, interacting with the M81 group. Prominent central emission features connected to a 
liner type activity prevent us from deriving the central kinematics (r < ±2 arcseconds) reliably. This galaxy exhibits 
an interesting shape of its velocity dispersion profile. The profile first rises gradually until « 150 kms~^ at radius 
of about 25 arcseconds. It then drops quickly to a minimum of about 130 kms^^ at 18 arcseconds and rises then 
again to « 160 kms^^ at the centre. The drop around 18 arcseconds is accompanied by a rapid change of slope of 
the rotational velocity which stays relatively flat outwards of this radius and a strengthening of the /13 moments. Also 
the otherwise vanishing moments rise to positive values (w 0.1) at r « 15 arcseconds where they form the most 
prominent double peak feature of our sample. The minor axis profile shows similar local minima in the dispersion 
profile that are accompanied by local maxima in the moments at radii of about 9 arcseconds. Given the inclination 
of 59° this points to a flattened structure within the bulge. 

NGC 3245 — .LARO*. — The HST F547M image reveals a nuclear dust disk while the larg er scale bulge morphology 



is classical. This galaxy potentially hosts a lens (de Vaucouleurs & de Vaucouleurs 1964). The rotational velocity 
exhibits a shelf reaching from about 3 arcseconds out to the bulge radius at 15 arcseconds. The velocity dispersion 
profile rises steadily towards the centre. The moments are anti-correlated with velocity inwards of 5 arcseconds. 
At 5 arcseconds they reach a maximum and then drop to zero towards the bulge radius. The /14 are compatible with 
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zero at the bulge radius (r^ = 8.5 arcseconds) but with decreasing radius briefly drop to about 
Finally inside of 2. arcseconds they become compatible with zero. 



.05 at ±5 arcseconds. 



NGC3521— .SXT4.. — The latest Hubble type with classical bulge in our sample. In HST F606W the classical bulge 
morphology stands in strong contrast to the strong outer disk spiral structure with a sudden transition of those two 
morphologies at r fa 10 arcseconds. The rotation curve forms a shelf at about 3 arcseconds but rises again slightly 
towards larger radii and starts to flatten out at 16 arcseconds. The velocity dispersion reaches a local maximum of 
K, 115 kms^^ at the bulge radius (r?, — 10.8 arcseconds) but then strongly drops with decreasing radius and reaches 
a local minimum of ~ 100 kms~^ at 4 arcseconds. Further inwards the dispersion is centrally peaked. Outside the 
bulge region the dispersion shows another local maximum around 40 arcseconds which also corresponds to a slight 
secondary shelf on the rotational velocity. The moments arc anti-correlated with velocity and form a shelf at about 
the same radii where the inner shelf in rotational velocity is seen. Outside the bulge /i3 becomes rather strong (up to 
0.22) which is accompanied by strong, positive /14 moments of up to 0.15. These l arge moments a r e a co nsequence of 
the double-peak structure of the LOS VDs which has been previously reported by Zeilinger et al. (2001 1 caused by a 
secondary kinematic component (see § 5.7). The minor axis velocity dispersion prohle does not show the strong local 
minima that the major axis profile shows but is similarly disturbed outside the bulge. 

NGC3898 — .SAS2.. — Outer disk dust spiral transitions into weak dust lanes over a smooth and classical bulge 
morphology in HST F606W. The rotation curve rises slowly out to r sa 10 arcseconds and flattens at about the 
bulge radius of « 15.7 arcseconds. The velocity dispersion profile is strongly centrally peaked and rises smoothly 
from sa 140 kms~^ to « 220 kms~^ in the centre. The ^13 moments are anti-correlated with velocity and reach local 
maximum at r « 10 arcseconds, the same radius where the flattening of the rotation curve sets in. 

NGC3992 — .SBT4.. — This galaxy is prominently barred. While the bulge morphology generally appears classical 
it exhibits a few randomly distributed dust lanes in the only available optical HST image in F547M. With increasing 
radius the velocity profile reaches a first plateau at r « 7 arcseconds but outside the bulge radius of 13.2 arcseconds 
it starts rising again slowly. The velocity dispersion profile is significantly depressed within the bulge region. The /13 
moments are strongly anti-correlated with velocity inwards of 7 arcseconds but become weaker further out. The /14 
moments show double peak feature around 8 arcseconds. 

NGC4203 — .LX.-*. — This galaxy represents a borderline case in the morphological bulge classification. T he bulge 
has a few dust lanes in HST F555W superimposed on a generally smooth morphology. Fisher & Drory (20081 classify 
the bulge as classical. The major axis velocity profile flattens out at about the bulge radius of 14.7 arcseconds. The 
velocity dispersion profile shows a prominent rise from about 105 kms~^ at the bulge radius to 175 kms~^ in the 
centre. The /13 moments are mostly compatible with zero in the bulge region. The /14 moments are noisy and show no 
significant trend within the covered region. 



NGC426O — .SBSl.. — Strongly barred. The bulge morphology is generally smooth. Martini et al. ( 2003 1 acknowledge 
the presence of dust structures in t he central region but c laim that they do not imply a sense of rotation. The bulge 
is classified as a classical bulge by jFisher fc Drory ( 2008 ) . The velocity dispersion profile is somewhat irregular and 
S-shaped within the bulge region, i'he /13 moments are small and mostly compatible with zero, The moments show 
a gradual drop from a central value of zero to about —0.06 at about 10 arcseconds. 

NGG4379 — .L..-P* — Exhibits a very smooth and featureless, classical central morphology in HST F555W. The 
bulge radius is 8.6 arcseconds. The inset of the flattening of the rotational velocity occurs already at 8-10 arcseconds. 
The velocity dispersion profile is centrally peaked, rising from w 90 kms~^ at 10 arcseconds to ss 120 kms""'^ in the 
centre. The /13 moments are generally small (within the errors mostly compatible with zero) but an overall trend 
points to anti-correlation with velocity. The /14 values are are noisy but mostly compatible with zero. 

NGC4698 — .SAS2.. — HST F606W shows some weak d ust lanes in the c e ntral re gion w hich ho wever do not i mply 
any sense of rotation. It is classified as classical bulge by Fisher & Drory (2010). |Falc6n-Barroso et al 
that the stellar velocity field displays rotation perpendicular to the major axis within the central 



( |2006[ ) find 
±5 arcseconds. 

The major axis rotational velocity is very slowly rising indicative of counter rotation. The velocity dispersion profile 
is mostly flat within the bulge region (r^ = 10.7 arcseconds) with a weak central peak. In the small radial range that 
our data cover the ft.3 and /i4 moments are mostly compatible with zero. 

NGC 4^72 — .SASl.. — - The bulge mor phology is smooth and featureless i n HST F6 06W, the bulge is consequently 
classified as classical by Fisher & Drory ( |2008 ). The [Oiii] map of ^ Falcon-Barroso et al. (2006|) s hows almost counter 



±5 arcseconds (see also Haynes et al. 2000) 



Our data virtually 

The velocity dispersion profile is generally noisy with values between 100 kms~^ and 



rotation with respect to the stars in the central 
only cover the bulge region. 

150 kms^^ but features a clear central peak. The /13 moments are anti-correlated with velocity and reach values of up 
to about ±0.07. /i4 moments scatter around values of 0.05 at all radii. 



Pseudohulges 

Possibly a triple barred system ( Erwin &: Sparke|1999 1. In HST F555W a dust spiral is seen, 



NGC2681— PSXTO.. 

which extends all the way into the centre. The rotational velocity curve shows a shelf between 



2 and sa 12 arcseconds. 
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The outer radius of the shelf region coincides with the bulge radius (r^ — 13.2 arcseconds). Inwards of 2 arcseconds 
the rotation curve drops quickly to zero. With decreasing radius, the velocity dispersion rises from about 50 kms^^ to 
a values of w 75 kms~^ at the bulge radius. Towards smaller radii it first stays relatively constant but then shows a 
step at about 5 arcseconds and rises again inwards of 2 arcseconds. The moments are anti-correlated with velocity 
in the region of the fast velocity rise but become correlated in the velocity shelf region. The /14 moments show a double 
peak feature in the radial range of 2-4 arcseconds. 

NGC2964— .SXR4*. — High contrast dust spiral in HST F606W. The bulge is small {n = 3.1 arcseconds) and 
we do not sufficiently resolve it to include this galaxy in any of our structural plots. Here we publish the kinematic 
profiles. The velocity dispersion rises from values of about 40 kms~^ in the disk to « 106 kms~^ at the bulge radius 
and exhibits a depression within the bulge, and /14 moments scatter strongly within values of ±0.1 which is possibly 
a consequence of the dust. 



NGC3166— .SXTO 
way into the centre 



. — The HST F547M im age shows strong dust features and a dust spiral which extends all the 
Laurikainen et al. ( 2004 1 describe this galaxy as strongly barred. The velocity dispersion profile 
shows a strong depression at « 4.5 arcseconds which is accompanied by relatively strong /13 moments of ±0.17 and 
positive hi moments and a local maximum in the rotation curve. The minor axis shows small but significant rotation 
within the bulge region (sa ± 10 kms^^) while the /13 moments are largely compatible with zero. The depression of 
the velocity dispersion is also seen on the minor axis whilst not as strong. 



NGC3351 — .SBR3.. — The bulge hosts a clear spiral structure and a nuclear ring f Fisher & Droryl 2010 1 and shows 



signs of active star formation. SAURON data ( ^Dumas et al.),2007j ) show a drop in the gas velocity dispersion derived 
from H/3 and lowered [Oiii]/H/3 ratios in the ring indicative of star formation. The velocity dispersion profile shows 
a depression within the bulge region — most strongly at r « 5 arcseconds down to a value of 70 kms~^ — and 
rises again centrally to « 90 kms^^ which is still below the values of about 100 kms~^ seen just inside of the bulge 
radius (r;, = 12.9 arcseconds). The moments are clearly anti-correlated with velocity. The minor axis profile shows 
significant rotation {vmax ~50 kms~^) indicative of a slit misalignment, in fact there is a 25°difference between our 
major axis slit position of 165°and the Hyperleda published value of 9.9°. 



NGC3368 — .SXT2.. — Complex mor phology with a number of stellar components ( Erwin|2004 Nowak et al.|2010 ) 



This galaxy is po ssibly double-barred (Jungwiert et al 



disk ( Erwin|2004 ) and is classified as pseudobulge by I' lsher &; Urory (2010). The bulge 
complex morphology is reflected in the kinematic structure, 'i'he rotational velocity reac 



iher & . 



The bulge hosts a strong nuclear spiral and an inner 



radius is 20.4 arcseconds. The 

iplex morphology is reflected in the kinematic structure, 'i'he rotational velocity reaches a local maximum at about 
7 arcseconds. This is accompanied by local minimum in the velocity dispersion profile that has been rising inwards 
until w 13 arcseconds. Inwards of 7 arcseconds the dispersion rises again, but asymmetrically about the centre, see 
also ^Nowak et al. (20101. The local maximum in velocity and local minimum in dispersion coincide with strengthened 
/13 moments. 7i3 is anti-correlated with velocity inwards of 15 arcseconds but correlated outside. The moments 
are close to zero at 13 arcseconds but become positive and reach a local maximum at about the same radii where the 
local maxima in velocity are observed and the drops in velocity dispersion and strengthening of /13 moments occur. 
The minor axis profile shows similar depressions in velocity dispersion. The minor axis velocity profile shows a central 
peak of about 30 kms~^. Visual inspection of the HET pre-acquisition images reveals that a minor offset of the slit 
position 1 arcsec) to the west is responsible for the peak. 

NGC3384 — .L BS-*. — Contains a nuc lear bar ([Fis her & Drory"2010') and a rapidly rotation disk described by 



Busarello et al. (1996); Fisher (1997) and 



Emsellem et al., (,2004) already find strong /13 moments in anti-correlation 



with velocity. 1 he velocity dispersion profile changes slope at 10 arcseconds and becomes more shallow towards smaller 
radii, b ut then exhibits a pron ounced peak inwards of 3 arcseconds that is accompanied by a dip (as seen by Fisher 
(1997); Emsellem et al. (2004) as well) in the /i4 moments that become positive just outside of this region. 



NGC3627 — .SXS3.. — Prominently barred galaxy with wide, open arms, inte racting with the Leo group. High 
contrast dust lanes that extend to the very centre are seen in HST F606W and let Fisher & Drory ( 2008 ) classify this 
as a pseudobulge. After a fast rise the rotational velocity forms a shelf between 3 arcseconds and the bulge radius 
of about 11 arcseconds. Towards larger radii the velocity rises again. The velocity dispersion rises inwards, starting 
already far outside the bulge radius. At ~ 4 arcseconds it flattens out and stays essentially constant. The /13 moments 
are anti-correlated with velocity inside of 9 arcseconds but change sign at larger radii and become correlated with 
velocity. While the minor axis rotation is compatible with zero at larger radii it exhibits significant rotation inwards 
of about 7 arcseconds that is also seen in anti-correlated ft.3 moments. 

NGC3675~ .SAS3.. — The HST F606W image clearly shows a high contrast flocculent spiral that extends aU the 
way into the centre. The velocity dispersion reaches a maximum of about 110 kms""'^ at « 40 arcseconds and stays 
rather constant inside of this radius. The /13 moments are anti-correlated with velocity within the bulge and reach 
maximum values of about 0.1 at the bulge radius of 8.5 arcseconds. 



NGC3945— RLBT+.. — Double barred ( |Kormendy||1979[ |1982[ [Wozniak et al.||1995[ [Erwin fc Sparke||1999 |2003, 

Erwin|2004 ) galaxy with prominent outer ring. Exhibits comjDhcated kinematic structure, 'i'he rotational velocity has 
local rninima around 18 arcseconds, rises then towards smaller radii and reaches a local maximum around 8 arcseconds 
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before it falls off to the centre. The dispersion profile has very strong local minima — drops from 150 kms ^ to 



110 kms" 



at r w 8 arcseconds but then rises again towards the centre. This galaxy shows exceptionally strong /13 



moments of up to ~ 0.25 which are anti-correlated with velocity in the inner region but become positively correlated 
at about 20 arcseconds. The /14 moments are similarly strong (up to 0.2) with significant central depression. The 
LOSVDs do show significant l ow v elocity shoulders at radii between 2" < ||r|| < 10" which is indicative of a kinematic 
distinct component (see Sec. 5.7). The minor axis profile also shows a local depression in velocity dispersion but 
at r w 3.5 arcseconds, the central dip in hi is seen as well. We measure slight rotation (i w ± 25 kms~^) along the 
observed position angle oi i = 64° and a significant central offset of the velocity (sa 40 kms^^). The major axis velocity 
profile shows that a slit position offset of one arcsecond is sufficient to explain the central velocity peak. A visual 
inspection of the pre-acquisition image confirmed that such an offset was indeed present (about 0.6 arcseconds). 

NGC4030— .SAS4.. — The flocculent spiral structu re — easily seen in H S T F606W — extends a ll the way to 
the centre with star-forming knots in the inner disk (Eskridge et al. 
pseudobulge 



2002). Fisher & Drory (2008) classify it as 



The major axis rotational velocity reaches «100 kms " at aDout ±iu arcsecona s. ine velocity aisp ei 
stays moderately flat within the bulge region and drops off outside , this was also observed by ("Ganda et al.|2 006) . The 

7i 



at about ±10 arcsecond s. The velocity disp ersion 
his was also observed by ("Ganda et al 
/13 moments are clearly anti-correlated with velocity and reach a values of up to 0.1 at a radius ot about 7 arcseconds. 
/14 moments are compatible with zero at all radii that are covered by our data. The bulge is small (rf, = 3.0 arcseconds) 
and we do not sufficiently resolve it to include this galaxy in any of our structural plots. 



NGC4274— RSBR2.. — Double barred galaxy (Sha w et al.||1995[ [Erwin||2004[ ). It is classified as pseudobulge with a 
Sersic index of 1.60 ± 0.35 by (Fisher &: Dr ory 2010|. 'I'he bulge hosts a prominent nuclear spiral including strong 
dust lanes and a nucle ar ring. The ring is se en as a fast-rotating, low-dispersion component in SAURON velocity and 
dispersion maps (Falcon-Barroso et al.|2006 ) . It also appears in their ionized gas maps through increased II/3 emission 
and lowered [OiiiJ/H/y ratios which indicates star-formation. The major axis rotational velocity rises quickly with 
increasing radius and starts flattening out at about 4 arcseconds, well within the bulge radius of 11.3 arcseconds. The 
velocity dispersion profile shows a strong depression inside the bulge with values of about 100 kms^^. Outside the 
bulge the dispersion rises to values exceeding 130 kms^^. The moments are well anti-correlated with velocity in 
the radial range that is covered by our data. They reach values of up to ±0.14 at about the same radius where the 
rotational velocity starts to fiatten out. This is accompanied by peaks in the moments with values of up to 0.14. 

NGC4314 — .SBTl.. — Strongly barre d galaxy with prominent nuclear ring. The ring appears in SAURON ionized 



gas maps (Falcon-Barroso et al. 2006) through lowered velocity dispersions, increase d H/3 and lowered [Oi iij/H/? 
which indicates star- format ion. 'Ihe bulge is consequently classified as pseudobulge by Fisher fc Drory| ( 2008 ) . The 
rotational velocity starts to flatten at the bulge radius of 8.6 arcseconds. The velocity dispersion profile is asymmetric 
but relatively flat, it varies between values of 105 kms^^ and 130 kms^^. The S/N only allows us to derive and hi 
moments inside a radius of 6 arcseconds. The /13 moments are mostly compatible with zero, the moments rise from 
zero in the centre to values around 0.05 at 4 arcseconds. 



Erwin & Sparke (1999) find a bright stellar ring that is notable 



NGC437I— .LBR+.. — Strongly barred galaxy.^ 

by adjusting the contrast of the e.g. the HST F6d6W image carefully or through unsharp masking. While free of 
obvious dust or s piral structures the ring with a radius of about 5 arcseconds falls within the bulge radius which lets 
Fisher & Drory ( |2010 ) classify this as a pseudobulge. The major axis rotational velocity starts to flatten out at a 
radius of 8 arcseconds — well inside the bulge radius or 22.9 arcseconds — with a weak shelf around 3-8 arcseconds. 
The velocity dispersion rises from about 105 kms~^ at the bulge radius to about 130 kms""'^ at 7 arcseconds. Inside 
of 7 arcseconds the velocity dispersion stays relatively constant. The ^,3 moments are somewhat asymmetric but 
generally anti-correlated with velocity within the bulge. The h^ moments are compatible with zero in the centre but 
rise gradually to values of about 0.05 at 10 arcseconds. The minor axis velocity dispersion stays mostly constant within 
the covered radial range. The minor axis /13 moments remain compatible with zero while the h^ moments rise from 
zero at the centre to values of about 0.05 at ±9 arcseconds. 



NGC4394 — RSBR3.. — Strongly barred galaxy with a face-on spiral in the central r w 5 arcseconds. |Fisher fc Drory| 
m08) classify it as a pseudobulge. The major axis velocity profile rises quickly to about 50 kms ^ at 3 arcseconds. 
'i'he velocity dispersion exhibits two prominent maxima around ±7 arcseconds. The maxima reach values of about 
105 kms~^ but drop quickly to about 80 kms~^ at 4 arcseconds and stay relatively constant towards smaller radii 
from there on. The /13 moments are anti-correlated with velocity inside the bulge. The /14 values are asymmetric and 
somewhat large within the bulge with values of up to 0.07 but we note that with dispersion values around 80 kms~^ 
we reach the limit of our ability to resolve those values properly. 



NGC4448-- .SBR2 
it as a pseudobulge 



Fisher & Drory (2010) note a mild spiral structure that extends into the centre and classify 
The major axis rotational velocity starts to flatten out at about 5 arcseconds — just inside of 
the bulge radius of 8.5 arcseconds. Outside of 5 arcseconds the velocity proflle exhibits a weak shelf. The velocity 
dispersion is mostly constant at a value of ~ 115 kms~^ at all radii covered by our data. The /13 moments are somewhat 
asymmetric, close to zero at positive radii but vary strongly at negative radii. The /14 moments scatter about zero in 
the bulge region with a few outliers at -0.05. 
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NGC4501— .SAT3.. — A nuclear spiral extends all the way to the centre in HST F606W. [Fisher fc Droryj (12010 1 
classify the bulge as pseudobulge. The major axis velocity profile rises quickly from the centre to the bulge radius of 
6.2 arcseconds but then flattens out at the bulge radius and forms a shelf out to about 20 arcseconds. The velocity 
dispersion rises from 75 kms""'^ in the disk to about 150 kms^^ at 15 arcseconds. The /13 moments are strongly 
anti-correlatcd with velocity inside the bulge and reach values of ± 0.1 at the bulge radius. /i4 moments are mostly 
compatible with zero at all radii. The minor axis velocity dispersion profile starts rising with decreasing radius inwards 
of 30 arcseconds and reaches a maximum of « 155 kms~^ at r = 10 arcseconds. The minor axis velocity dispersion 
exhibits a central depression of about 15 kms^^. The minor axis velocity, /13 and /14 moments are mostly compatible 
with zero. 



NGC4536 — .SXT4.. — A strong dust spiral extends into the very centre. Fisher & Drory (2008) classify the bulge 
as pseudobulge. The bulge radius is 10.1 arcseconds. The major axis rotational velocity jfiattens well inside bulge 
radius at around 2 arcseconds. The velocity dispersion profile is slightly asymmetric but mostly constant within the 
bulge region. The moments are anti-correlated with velocity inside the bulge region and reach values of about 0.1. 
The /14 profile is asymmetric. The minor axis velocities show slight asymmetric rotation (ss 25 kms~^). The velocity 
dispersion profile rises centrally but stays relatively flat within ± 10 arcseconds. The minor axis moments are 
somewhat noisy but seem to show a central depression. 

NGC4569 — .SXT2.. — A nuclear spiral extends all the way to the centre. The bulge is classified as pseudobulge by 



Fisher & Drory (2010). The major axis rotational velocity rises with increasing radius to a local maximum of about 
80 kms~^ at ± 3 arcseconds. The bulge radius is — 9.6 arcseconds. The velocity then drops to about 50 kms^^ at 
the maximum radius covered by our data. The velocity dispersion rises with decreasing radius to about 100 kms^^ 
at a radius of 6 arcseconds and then drops and reaches a local minimum at around 3 arcseconds, roughly coinciding 
with the locations of the local maxima in the velocity profile. The moments are anti-correlated with velocity inside 
the bulge. Their absolute values reach up to 0.15. The /14 moments show a strong double peak feature at about 
3 arcseconds and fall off to zero at the bulge radius. The minor axis profile shows rotation in the bulge region. While 
somewhat asymmetric, the minor axis dispersion profile does not show the same complicated structure of the major 
axis profile. The moments on the minor axis are mostly compatible with zero. The ^13 moments on the minor 
axis are mostly slightly positive with a mean value of 0.03. The /14 moments are generally noisy but the double peak 
feature of the major axis is reproduced. 



NGC4736 — RSAR2.. — Hosts a nuclear bar ( Sakamoto et al.||1999 ) and prom inent nuclear sp iral which extends all 



the way into the centre in HST F555W. The bulge is classified as pseudobulge by Fisher & Drory (2010). The obtained 
kinematic data extend well into the disk. The rotational velocity flattens out abruptly at about the bulge radius of 
14.2 arcseconds and shows a shallow negative gradient out to about 70 arcseconds where our data points start to 
become sparse. The velocity dispersion rises abruptly from about 75 kms^^ to 115 kms~^ at about the bulge radius. 
Well within the disk at radii larger than 50 arcseconds we see again a gradual increase of velocity dispersion. Inside of 
2.5 arcseconds the velocity dispersion exhibits a central drop. The moments are anti-correlated with velocity but 
show s-shape around the centre. They reach exceptionally large values of ±0.2 at the bulge radius. The /14 moments 
are compatible with zero in the inner bulge but reach pronounced local maxima of values as large as 0.25 at about 
the bulge radius. They fall off to zero at r sa 35 arcseconds. These strong highe r moments are a consequence of the 



multi-component structure of the LOSVDs at the respective radii (see Sec. 5.7 1. The minor axis profile reflects the 
rich structure seen in the major axis profile. The velocity dispersion rises significantly inwards of 10 arcseconds. The 
/13 moments are mostly compatible with zero at all radii, moments are zero inside of 10 arcseconds but rise to about 
0.1 at 20 arcseconds. 

NGC5055 — .SAT4.. — The HST F606W image shows that the outer disk flocculent spiral extends into the very 



centre. Consequently this galaxy is classified as pseudobulge by Fisher & Drory (2010). The velocity dispersion profile 
remains fiat inside the bulge radius of 18.3 arcseconds. The /13 moments are anti-correlated with velocity inside the 
bulge and reach absolute values of up to 0.1 at the bulge radius. The moments are compatible with zero inwards 
of 10 arcseconds, they become noisy further out but show a weak tendency towards more positive values towards 
the bulge radius. The minor axis velocities appear somewhat irregular but small (< 20 kms~^). The minor axis /13 
moments are noisy but mostly scatter close to zero. Again the /14 moments are mostly compatible with zero inwards 
10 arcseconds. but show a weak increase further out but only on the east side. 

NGC5248~ .SXT4.. — Has a prominent nuclear spiral clearly visible in HST F814W. SAURON maps show the 
presence of a nuclear ring in H/3 and [O ill] emission. A l owered [O iii]/H/3 shows that the ring is star-forming. The bulge 
was classified as pseudobulge by Fisher & Drory (2 010| ) . The rotational velocity starts to fiatten at about 6 arcseconds 
— well inside the bulge radius of 15.4 arcseconds. 'ihe velocity profile shows a shelf between sa 10 arcseconds and 
w 40 arcseconds. The velocity dispersion profile is mostly flat at about 80 kms~^ with two small peaks at r w 
±4 arcseconds. The moments are anti-correlated with velocity and reach values of up to ±0.1. The /i4 moments 
scatter around values of 0.05. 

NGC5566 — .SBR2.. — Shows a nuclear spiral in HST F606W. The surface brightness profile does not resemble a 
traditional bulge plus disk structure. We do not include this galaxy in any of the structural plots and publish only 
the kinematic profile here. The rotational velocity starts to fiatten at 6 arcseconds. The velocity dispersion profile 
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rises towards the centre and peaks at a value of about 150 kms^^. The moments are anti-correlated with velocity 
and reach values of up to 0.15. The /14 moments are small in the central arcseconds but rise to values of about 0.1 at 
~ ±5 arcseconds. 



N GCini — -SXRB . . — A nuclear bar extends out to about r = 10 arcseconds. The bulge is classified as pseudobulge 



by Fisher & Drory (2010). The flattening of the major axis rotational velocity coincides with bulge radius of = 
8.6 arcseconds. The velocity dispersion rises from about 50 kms^^ in the disk to values of w 115 kms^^ inside the 
bulge but remains relatively flat inside the bulge radius. The moments are weakly anti-correlated with velocity but 
remain small. The moments drop to values of -0.05 at « ±5 arcseconds. The minor axis shows an asymmetric 
velocity profile with values of up to ±25 kms~^. The minor axis velocity dispersion again rises from about 50 kms~^ in 
the disk to values of about 115 kms^^ at 2.5 arcseconds and remains flat inside. The /13 moments are anti-correlated 
with velocity. The /14 moments are noisy and scatter around zero. 



NGC 1H3 — RLBS-I-.. — The central region exhibits some weak dust la nes overlaid on a gener ally (Martini et al. 
20031 smooth light distribution. The bulge is classified as classical bulge by Fisher & Drory (2008). The amplitude of 
the rotation is small at about ±25 kms~^ due to the low inclination. We find rotation of similar value along the minor 
axis slit due to a misplacement of the slit (166° rather than 10°). The velocity dispersion is flat for both position 
angles and takes values of about 80 kms~^. Both, major axis and minor axis /is moments are anti-correlated with 
velocity and become compatible with zero at the bulge radius (r?, — 5.6 arcseconds). The /14 moments are very noisy, 
probably due to the low velocity dispersion of this object. 



Bulges without classification 

NGC246O — .SASl.. — Mixed type morphology in HST F606W with a weak asymmetric dust structure in the bulge 
region that is overlaid on an otherwise smooth light distribution. We label this bulge as unclassified. The decomposition 
gave a value of 3.5 ± 0.32 for the Sersic index and 6.6 arcseconds for the bulge radius. Within this region the velocity 
dispersion profile is flat and /13 and /14 moments scatter around zero. 



NGC3593 — .SASO*. — Peculiar bulge s t ructu re with prominent spiral visible even in NIC F160W. The bulge is 
classified a pseudobulge by Fisher & Drory (2010). We label it as non- classified because its high inclination inhibits an 
unperturbed view into the bulge region. Ihis is the only example in our sample where counter rotation is seen in the 
vel ocity profile as an ac tual change of the sign of the velocity with respect to the systemic velocity (this was found also 



by Bertola et al. 1996 ) . The counter rotation within the bulge radius is reflected in the anti-correlated ft.3 moments. 
The velocity dispersion drops dramatically from « 115 kms^^ at the bulge radius to « 60 kms^^ in the centre. 



NGC3953 — .SBR4.. — Fisher & Drory (2010) classify this bulge as pseudobulge but acknowledge that there is no 
optical HST image available. Here we label it as non- classified. The rotational velocity profile first reaches a weak 
local maximum at r « 2.5 arcseconds before it starts rising again outside of 8 arcseconds. The disk velocity dispersion 
rises centrally from values of « 50 kms~^ at ± 30 arcseconds to « 130 kms^^ at « 8 arcseconds. Inside of a radius of 
6 arcseconds it then falls toward a central value of « 110 kms~^. The fast central increase of velocity is accompanied 
by strong anti-correlated /13 moments with values of up to ± 0.1 at r 6 arcseconds, they become correlated with 
velocity outside of 8 arcseconds. The moments are generally noisy in the bulge region and scatter between zero and 
0.05. 



NGC4826 — RSAT2.. — Also named the hlack eye galaxy. An extreme dust spiral in the central 50 arcseconds stands 
in strong contrast to a vi rtually dust free outer disk. The central dust content leads to a classification as pseudobulge 



Fisher & Drory| (|2008|). The major axis kinematic profile is rich in structure. 
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The rotational velocity rises quickly 

from the centre to a value of ~ 50 kms~^ at r = 4 arcseconds. It then forms a shallow trough around 8 arcseconds and 
then rises again — more slowly — out to 50 arcseconds where it finally flattens out. The velocity dispersion in the disk 
is « 45 kms""'^, it shows a distinct central increase inwards of 50 arcseconds. The dispersion reaches values of up to 
110 kms~^ inside of the bulge radius (r;, = 25.4 arcseconds). From there on it stays relatively constant with decreasing 
radius except for a mild depression down to 90 kms~^ in the central few arcseconds. The /13 moments are strongly 
anti-correlated with velocity for r < 9 arcseconds and reach absolute values of up to 0.15. The /i4 moments show two 
peaks at about ±3 arcseconds. We find weak rotation along the minor axis (~ ±10 kms~^). The central increase in 
velocity dispersion is also seen along the minor axis, the increase sets in at a radius of about 25 arcseconds. This is 
much closer to the centre than in the case of the major axis and points to a flattened structure. For the photometry we 
obtain a mean bulge ellipticity of 0.23. The radial difference of the dispersion increase along the major and the minor 
axis would point to an ellipticity of about 0.45. The centre of the velocity dispersion profile is asymmetric which may 
be a consequence of the strong dust. The minor axis /13 and moments are compatible with zero. The fact that the 
final flattening of the rotational velocity and the inset of the dispersion increase at 50 arcseconds falls together with 
the sudden appearance of the strong dust structure is intriguing and lets one suspect that the actual bulge radius fall 
closer to 50 arcseconds. The much smaller bulge radius fr om t he decomposition may be a consequence of the strong 
dust content in the central regions (see also discussion in § 5.4 1. 
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NGC72n — RSAR2.. — Fisher & Drory (2010) point out a sudden break in morphology at a radius of about 
8 arcseconds where the outer spiral transitions onto a relatively smooth morphology with little dust. They consequently 
label this bulge as classical. However we determine a bulge radius of 11.2 arcseconds where there is already a pronounced 
spiral pattern visible. We label this galaxy as non- classified. The rotational velocity profile starts to flatten out at 
about the bulge radius of. The major axis dispersion profile is asymmetric with higher values on the east side of the 
centre. Within the central ± 2 arcseconds the velocity dispersion shows a mild depression of about 20 kms~^. The 
/i3 moments are well anti-correlated with velocity. The /14 moments scatter around values of 0.05. We observe mild 
rotation on the minor axis (of the order of ±10 kms~^) indicative of a slight slit misalignment. The minor axis velocity 
dispersion is symmetric and rises towards the centre from values of about 100 kms~^ at radii of ± 20 arcseconds to 
145 kms""'^ at 2 arcseconds. Within the central arcseconds the mild depression which is seen on the major axis is 
reproduced on the minor axis. The minor axis moments mostly scatte r around zero while ^4 mome nts fall closer to 
0.05 with a few relatively large outliers at radii around 6-10 arcseconds. Merrifield & Kuijken (19941 found the 20%- 



30 % of the light is captured in a counterrotating component. We confirm this and present a kinematic decomposition 
^5Jl 



NGC7331 — .SAS3.. — The HST F555W image shows several dust lanes in the bulge region. H owever these do 
not im ply any sense of rotation and seem to be overlaid on a generally smooth light distribution. Fisher fc Drory| 



(2008) admit that the high inclination leaves the classification as classical bulge questionable. Here we label it as 
non- classified. The rotational velocity profile is already flattened at the bulge radius of rf, — 26 arcseconds. The 
velocity dispersion rises from about 75 kms^^ in the disk to 125 kms~^ in the centre. The dispersion profile has two 
steps or shoulders at w ±20 arcseconds. The /13 moments are generally anti-correlated with velocity and reach a local 
maximum of about ±0.15 at r = 15 arcseconds. Also the /14 moments reach local maxima with values of up to 0.15 in 
the same radial range. These large moments are a cons equence of the doub le-pea k str ucture of the LOSVDs caused by 
a counterrotating kinematic component discovered by Prada et al. ( 1996[ ) (see § 5.7). 



